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An  analysis  of  numerical  methods  for  extracting 
aerodynamic  coefficients  from  dynamic  test  data  has  been 
conducted.   The  emphasis  of  the  analysis  is  on  the  effects 
that  random  measurement  errors  in  the  data  and  random 
disturbances  in  the  system  liave  on  the  accuracy  witli  which 
the  coefficients  for  linear  and  nonlinear  systems  can  be 
determined.   Both  deterministic  and  stochastic  methods  for 
extracting  the  coefficients  and  determining  their  uncertain- 
ties are  considered. 

The  deterministic  technique  considered,  due  to  Chapman 
and  Kirk,  provides  excellent  estimates  of  both  linear  and 
nonlinear  static  pitching  moment  coefficients  for  the  range 
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of  measurement  errors  and  system  noise  considered.   Someivhat 
less  accurate  estimates  of  linear  damping  coefficients  are 
obtained.   Nonlinear  pitch  damping  coefficients  extracted 
using  this  deterministic  technique  are  affected  considerably 
by  both  measurement  errors  and  system  noise.   The  estimated 
standard  deviations  of  the  extracted  coefficients  obtained 
using  standard  techniques  are  generally  adequate  when  the 
data  being  analyzed  contain  only  measurement  errors. 

The  stochastic  approach  considered  demonstrates  the 
feasibility  of  using  an  extended  Kalman  filter,  with  a 
parameter  augm.ented  state  vector,  for  determining  the  values 
of  the  aerodynamic  coefficients  and  their  uncertainties 
from  noisy  dynamic  test  data. 

The  specific  filter  used  generally  reaches  near  steady- 
state  conditions  in  its  estimates  of  the  parameters  in  less 
than  one  second.   Variations  in  the  initial  parameter 
variances  or  in  the  estimates  of  the  noise  statistics 
essentially  affect  only  the  determination  of  the  nonlinear 
damping  parameter. 

Parameter  estimates  obtained  with  the  extended  filter 
compare  favorably  with  previously  obtained  results  using 
deterministic  techniques.   Estimates  of  the  parameter 
uncei tainties  provided  by  the  filter  are  generally  superior 
to  those  obtained  with  deterministic  techniques  particularly 
wlien  system  noise  has  corrupted  the  data. 
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CHAPTER  I 
INTRODUCTION 

This  dissertation  presents  an  analysis  of  methods 
for  determiming  aerodynamic  coefficients  of  flight  vehicles 
from  dynamic  test  data.   This  determination  is  accomplished 
by  finding  niimerical  values  of  the  aerodynamic  coefficients 
appearing  in  the  equations  of  motion  such  that  the  solutions 
to  these  equations  are  adequate  representations  of  the 
given  test  data.   When  this  determination  has  been  made, 
the  coefficients  are  said  to  have  been  extracted  from  the 
dynamic  data. 

The  emphasis  of  the  analysis  is  placed  on  the  effect 
that  random  iseasurement  errors  in  the  data  and  random 
disturbances  in  the  system  have  on  the  accuracy  with  which 
the  coefficients  for  linear  and  nonlinear  systems  can  be 
determined.   Both  deterministic  and  stochastic  system 
models  for  extracting  the  coefficients  and  determining 
their  uncertainties  are  considered. 

Historical  Background 

The  determination  of  aerodynamic  coefficients  from 
dynamic  data  is  generally  agreed  to  have  begun  with  the 


work  of  Fowler,  Gallop,  Lock  and  Richmond   in  the  first 
quarter  of  this  century.   Their  basic  technique  was  con- 
cerned with  determining  moment  and  damping  characteristics 
of  artillery  shells  by  firing  these  projectiles  through 
spaced  cards  and  reconstructing  the  pitch  and  yaw  angle 
time  histories  by  observing  the  obliqueness  of  the  holes 
that  resulted  when  the,  shells  passed  through  the  cards. 

Their  technique  of  data  measurement  is  still  in  use  at 

2 
some  ballistic  ranges  to  this  date. 

3 
Nielsen  and  Synge   later  clarified  the  linear  theory 

of  Fowler  et  al .  in  their  work  during,  and  immediately 
following,  Vvorld  War  II. 

Coefficient  extraction  techniques  that  are  currently 
in  use  at  many  ballistic  and  wind  tunnel  facilities  evolved 
from  the  work  of  Murphy  '   and  Nicolaides,   both  of  whom 
have  considered  various  combinations  of  degrees  of  freedom 
for  a  variety  of  flight  vehicles.   Their  aerodynamic  co- 
efficient extraction  techniques  feature  a  least  squares 
fit  to  dynamic  data  of  the  exact  solutions  of  linear 
equations  of  motion,  or  approximate  closed  form  solutions 

of  slightly  nonlinear  equations  of  motion  using  the  quasi- 

7 
linearization  technique  of  Kryloff  and  Bogoliuboff.    Ex- 


tensions of  the  work  of  Murphy  and  Nicolaides,  particularly 
for  more  coir.plex  nonii 
Eikenberry  and  Ingram, 


for  more  coir.plex  nonlinear  systems,  have  been  made  by 

9 


The  requirement  for  closed  form  solutions  of  the 
equations  of  motion  has  recently  been  eliminated  by  the 
formulation  of  least  squares  techniques  which  fit  numerical 
solutions  of  the  equations  of  motion  to  dynamic  data.   The 
minimization  of  the  least  squares  criterion  function  involves 
differential  corrections,  as  do  the  techniques  of  Murphy 
and  Nicolaides;  the  required  partial  derivatives  are 
determined  by  numerically  integrating  parametric  differ- 
ential equations  which  are  derived  from  the  equations  of 
motion.   Contributions  to  this  numerical  coefficient 
extraction  technique  have  been  made  independently  by 

Chapman  and  Kirk,    Knadler,    Goodman   '    and 

14 
Meissinger,  "  although  the  method  is  usually  referred  to 

as  the  Chapiaan-Kirk  technique.   The  computational  require- 
ments of  this  technique  are  sometimes  extensive,    but 
this  is  usually  outweighed  by  the  fact  that  it  can  be 

used  to  analyze  highly  nonlinear  aerodynamic  forces  and 

.      16 
moments . 

All  of  the  coefficient  extraction  techniques  that 
have  been  discussed  to  this  point  are  deterministic  in 
nature  in  that  the  modeling  of  the  equations  of  motion 
does  not  account  for  random  disturbances  in  the  system. 
Most  angular  or  translational  motion  data  obtained  from 
ballistic  raages  or  wind  tunnel  dynamic  tests  have,  never- 
theless, beea  affected  by  random  system  noise  and  random 


measurement  errors.   The  effects  of  noise  of  these  types 

on  the  accuracy  of  coefficients  extracted  from  dynamic  data 

17 
using  deterninistic  techniques  are  of  current  interest. 

9 
With  the  exception  of  Ingram's   partial  analysis  of  noise 

Q 

effects  on  Eikenberry's  "Wobble"  program,   however,  little 
has  been  reported  concerning  these  effects.   One  of  the 
primary  goals  of  the  research  reported  here  is  to  show  some 
effects  of  random  measurement  errors  and  system  noise  on 
the  accuracy  of  aerodynamic  coefficients  extracted  from 
dynamic  data  using  the  most  general  of  the  various  determi- 
nistic techiiiques,  that  of  Chapman  and  Kirk. 

In  recent  years  increasing  attention  has  been  given 
to  parameter  and  state  variable  estimation  through  the 
use  of  stochastic  modeling  of  the  physical  system  of 
interest.   'lOst  of  this  work  has  been  done  by  optimal 
control  specialists  and  is  a  direct  result  of  pioneering 

contributions  in  linear  filtering  theory  of  Kalman  ^' 

1  9 
and  Bucy.^    The  Kalman  filter  provides  estimates  of  the 

states  of  noisy,  linear  dynamic  systems  as  well  as  estimates 

20 
of  the  state  variable  uncertainties.   Bryson  and  Ho, 

among  others,  have  extended  the  Kalman  filter  theory  to 

include  estimates  of  the  states  for  nonlinear  systems. 

7 1 
Mehra'""   has  recently  proposed  a  maximum  likelihood 

technique  for  the  determination  of  aerodynamic  coefficients 

from  dynamic  data  using  the  Kalman  filter  for  linear  systems 


and  the  extended  Kalraan  filter  when  the  system  is  non- 
linear.  Results  verifying  his  proposals  are,  as  yet, 
unavailable. 

An  additional  maximum  likelihood  technique  has  been 

22 
developed  by  Grove  et  al .    and  has  been  used  with  modest 

23 
success  by  Suit. 

A  stochastic  approach  to  the  coefficient  extraction 

problem  using  an  extended  Kalman  filter  with  parameter 

augmented  state  vector  is  developed  in  this  dissertation. 

Results  obtained  with  this  filter  are  also  presented.  . 

Scope 

The  following  chapters  are  devoted  to  an  analysis 
of  the  effecti  of  random  system  noise  and  measurement 
errors  on  coefficients  extracted  from  dynamic  data  using 
the  method  of  Chapman  and  Kirk,  as  well  as  the  development 
and  evaluation  of  an  extended  Kalman  filter  for  solving 
essentially  the  same  problem.   The  equation  of  motion 
used  in  the  analysis  is  one  describing  the  one-degree- 
of-freedom  pitching  motion  of  a  rigid  body.   Pitching 
moments  that  are  both  linear  and  nonlinear  functions  of 
angle-of -attack  are  considered. 

The  theory  of  the  Chapman- Kirk  technique  is  recounted 
briefly  in  Chapter  II,  followed  by  tlie  analysis  of  noise 
effects  on  this  method  in  Chapter  III. 


Developisaents  leading  to  the  formulation  of  the  extended 
Kalman  filter  to  be  used  are  given  in  Chapter  IV.   Results 
obtained  with  the  extended  filter  are  presented  and  discussed 
in  Chapter  V. 

Some  concluding  remarks  are  offered  in  Chapter  VI. 

Computer  program  listings  and  definitions  of  program 
variables  are  given  in  the  two  appendices  following  the 
body  of  the  paper. 
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CHAPTER  II 


THE  COEFFICIENT  EXTRACTION  TECHNIQUE 
OF  CHAPMAN  AND  KIRK 


A  brief  reconstruction  of  tlie  Chapman-Kirk  algorithm 
for  extracting  aerodynamic  coefficients  from  dynamic  data 
is  given  in  this  chapter.   The  standard  technique  for 
estimating  uncertainties  in  coefficients  determined  from 
a  least  squares  fit  of  a  given  function  to  experimental 
data  is  also  presented.  ■ 

Development  of  the  Extraction  /Algorithm 

The  extraction  technique  of  Chapman  and  Kirk  has  two 
very  basic  requirements:  (1)  the  different  al  equation 
of  motion  for  the  body  of  interest  must  be  given  and 
(2)  a  set  of  experimental  data  based  on  the  observed  motion 
of  the  body  must  be  available. 

As  an  example,  consider  the  nonlinear  equation  of 
pitching  motion  for  a  rigid  body 

a  +  (C4  +  C5a^)i  +  (C3  +  C5a2)a  =  0       (1) 

subject  to  the  initial  conditions 

a(0)  =  C^ 
"^  i(0)  =  C2 
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J 

where  a  indicates  the  derivative  of  a  with  respect  to  time. 
Suppose,  also^  that  the  time  history  of  the  pitch  angle 
as  recorded  during  some  experiment,  a  (t-),  i=l,2,...m, 
is  also  available. 

The  teclhnique  of  Chapman  and  Kirk  is  used  to  determine 
the  values  o>£   the  C  (j  =1 , 2 ,  .  .  .  ,6)  in  Equation  (1)  which 
result  in  the  solution  to  this  equation  of  motion  being 
a  best  fit  to  the  test  data  in  a  least  squares  sense. 
Thus  it  is  necessary  to  minimize  the  least  squares  cri- 
terion function 


m 
i  =  l 


a3(t.)  -  a^(t.) 


(2) 


where  a  (t-/  is  obtained  from  the  solution  to  Equation  (1) , 
Now  it  is  well  known  that  in  order  to  determine  para- 
meters directly  by  a  least  squares  fit  of  a  given  function 
to  test  data,,  the  parameters  must  appear  linearly  in  the 
function.   Tiiis  requirement  is  met  in  the  problem  at  hand 
by  expanding  a      in  a  truncated  Taylor  series  about  the 
numerical  solution  resulting  from  some  initial  estimates 
of  the  parairc-sters  of  interest,  C-  ,  i.e.. 


-c^'-O    =    "co^^^i)  ^  .1 


3a 


AC 


(3) 


where 


da 


dC 


are  evaluated  for  C-=C.   and  AC.=C.-C.  , 

J   JO        J   J   JO 


Substituting  (3)  into  (2)  yields 


m 

S  =   I 
i  =  l 


e  ^  1"^ 


a   (t.) 
CO    1 


I 


' 

3a 

AC. 

J 

i        ^ 

(4) 


Now,  assuming  that  the. 


3a 


3C 


are  known,  Equation  (4)  is  a 


function  of  the  AC.'s  only.   Therefore,  to  determine  the 
values  of  these  coefficients  that  will  minimize  this  equation, 
it  is  necessary  to  take  the  partial  derivative  of  Equation 
(4)  with  respect  to  each  AC.,  set  each  of  the  resulting 
equations  to  zero  and  solve  for  the  AC.'s.   Carrying  out 
these  operations  yields  the  following  matrix  equation 


[A]  [AC]  =  [B] 


(5) 


where  A  is,  in  this  case,  a  6  X  6  matrix  with  elements 
given  by     . 


^Jk 


m 

=  I 
i-ll 


3a 


aCj 


3a 


3C 


k) 


C6) 


AC  is  a  6  X  I  column  matrix,  or  vector,  and  B  is  a  6  X  1 
column  matrix  with  elements  given  by 


B.  = 


m  r 

I 

i  =  l 


(t^) 


'^i' 


3a 
3Cj 


(7) 
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The  solutioji  to  Equation  (5)  is 

[AC]  =  [A]-l  [B]     .  (8) 


The  solinations  for  the  AC.'s  obtained  from  Equation  (8) 
are  exactly  correct  only  if  a   is  a  linear  function  of  the 
C-'s  as  assumed  by  Equation  (3).   This  condition  of  linearity 
is  seldom  the  case  and  the  process  must  be  repeated  with 
new  initial  guesses. 


C. ,  =  C.   +  AC.  (9) 

J  1      JO       JO  ^    ■' 


until  the  cikange  in  the  criterion  function  [Equation  (2)] 
from  one  iteration  to  the  next  is  sufficiently  small. 
The  algorithm  just  presented  requires  that  time 
histories  of   the  influence  coefficients,  ^7=r-,  be  available. 
These  time  hListories  are  determined  by  numerically  inte- 
grating para-Bietric  differential  equations  which  are  derived 
by  dif f erenc iating  the  equation  of  motion  with  respect  to 
each  of  the  parameters  of  interest.   As  an  example,  the 
parametric  differential  equation  for  Cj,  the  initial 
condition  of  a.    is 


■^[■d    +     (C4  +  C6a2)a    +     CC3+C5a2)ct]     =     0 


11 


-— -—  +  -— —  a  +   ^   a'^a  +  ZLrCta.    —77 

3Ci   aCi     aCi        s   3C 


6    3C- 


9C3 


3a     3C  5  3  a 

+  C3  -TTT-    + a3  +  3C5a2  —7-  =  0 

^3Li    3a  ^    3Li 


Assuming  that  the  parameters  are  independent  of  each  other 
and  that  the  order  of  differentiation  can  be  reversed, 
the  final  form  of  the  desired  equation  is 


df 


3a 


ISCiJ 


^e- 


dt 


3a 


UcJ 


+  (C3+3C^a2+2Cgaa) 


3a 


3C 


(10) 


subject  to  the  initial  conditions 


3a(0)  _  3C_L 
3C,   "  3C, 


=  1 


3a(0)     9C2 


3C 


3C 


=  0 


(11) 


The  complete  set  of  parametric  differential  equations  for 
the  given  equation  of  motion  [Equation  (1)]  is  given  in 
Chapter  III. 

In  summary,  the  process  for  extracting  numerical 
coefficients  from  test  data  given  the  system  model  [Equa- 
tion (1)]  and  criterion  function  [Equation  (2)]  is 

1.  Estimate  the  numerical  values  of  the  C.'s. 

2.  Integrate  Equation  (1)  to  obtain  c'co'-^i-'" 
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3.   Determine 


3a 


4.  Solve  Equation  (8)  for  the  AC.'s. 

5.  Repeat  the  process  with  C.  =C.  +AC.   until  the 

^        ^  J  ^   JO    JO 

chsBge  in  Equation  (2)  is  sufficiently  small. 

Estimation  of  Extracted  Parameter  Uncertainties 

The  estimation  of  the  uncertainties,  or  standard 
deviations,  of  parameters  that  have  been  determined  by  the 
least  squares  fitting  of  a  given  function  to  test  data  is 
a  well-known  result  available  in  a  variety  of  references 
(see,  for  exa^iiple.  References  24  or  25). 

The  leasr  square  parameter  uncertainties  are  estimated 
in  this  paper  by 


a.    =    /A .  .  -  1  .  /-^ 
J      JJ       m-K 


(12) 


w]iere  A..'^  is  the  j  th  diagonal  element  of  the  inverse 
J  J 

Grammian  matrix  [Equations  (6)  and  (8)],  S  is  the  sum  of 
the  squares  of  the  residuals  as  given  by  Equation  (2) ,  m 
is  the  total  number  of  data  points  and  K  is  the  total  number 
of  parameters  being  determined  by  the  fit. 


CHAPTER  III 

ANALYSIS  OF  THE  CHAPMAN- KIRK 
COEFFICIENT  EXTRACTION  TECHNIQUE 


Coefficient  Extraction  Computer  Program 

This  section  provides  a  detailed  description  of  a 
one-degree-of -freedom  coefficient  extraction  computer 
program  based  on  the  previously  described  iterative  process 
of  Chapman  and  Kirk.   This  program  considers  the  pitching 
motion  of  a  symmetric  missile  about  a  fixed  point.   It 
is  used  to  determine  values  of  static  pitching  moment 
coefficient  derivatives,  pitch  damping  coefficients,  and 
a  trim  term  so  that  the  solution  to  the  appropriate  differ- 
ential equation  of  motion  is  a  best  fit  to  test  data  in  a 
least  squares  sense. 

The  program  was  developed  to  be  used  as  an  economical 
tool  in  the  evaluation  of  the  sensitivity  to  noise  and 
convergence  sensitivity  of  the  Chapman-Kirk  technique  when 
operating  in  its  least  complex  mode. 

Computational  Equations 
Equation  of  motion 

The  complete  equation  of  motion  that  is  Ubcd  in  the 
program  is 


13 


14 


a    +     CC    +C    a2)a    +     (C    +C    a2+c    ct'*)  a    +    C„    =    0 
H         6  3         5  7  8 


(13) 


with    initial    conditions 


accy)   =  Ci 


a(0)     =    C. 


(14) 


where 


CCmao)qAd 
I 


C,     = 


(Cmqo)qAd^ 
2VI 


(Cjnaz)^^^ 


C      = 

6 


CCniq2)qAd- 
2VI 


^7     = 


(Cma4)qAd 


(Cm6a)qAd 


Parametric  differential  equations 

The  eiglit  parametric  differential  equations  for  the 
given  equation  of  motion  [Equation  (13)]  and  initial  con- 
ditions [Equation  (14)]  are  of  the  form 


P.  +  AP-  +  BP.  =  F 
J      J      3     . 


1,2, 


(15) 


where 


^j 


9a 
3C. 


d  da 


d2 


da 


dt  dC- 


dt2  3C. 


A  =  (C,  +  C,a2) 
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and 

B  =  ^C^    +    SCgO^  +  2Cgai  +  SC^a*^)  . 

The  values  of  the  initial  conditions  P.(0)  and  P  (0)  as 
well  as  the  functional  forms  of  the  nonhomogeneous  term 
F-  are  given  in  Table  I  for  j=l,2,...,8. 


Table  I.   INITIAL  CONDITIONS  AND  NONHOMOGENEOUS 

TERMS  FOR  PARAMETRIC  DIFFERENTIAL  EQUATIONS 


j  Pi  CO)  ^.(0)  F. 

110  0 

2  0  1  0 

3  0  0  .  -a 

4  0  0  -a 

5  0  0  -a3 

6  0  0  -act^ 

7  0  0  -a5 

8  0  0-1 
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Program  Description 

The  program  is  written  in  Fortran  IV  for  use  primarily 
on  an  IBM  360/65.   The  paragraphs  that  follow  describe  the 
functions  of  the  main  program  and  its  four  subroutines, 
the  required  input  data,  and  the  program  options.   A  listing 
of  the  complete  program  is  given  in  Appendix  I. 

Main  program  functions 

The  functions  of  the  main  program  in  their  approximate 
order  of  occurrence  are  as  follows: 

(a)   Reads  and  writes  input  data. 

Cb)   Computes  initial  parameter  values  (these  are 
the  C.'s  that  appear  in  Equation  (13). 

(c)  Calls  the  numerical  integration  subroutine  ADDUM. 

(d)  Writes  current  parameter  values. 

(e)  Computes  the  matrix  elements  required  for  incre- 
menting the  parameter  values. 

(f)  Calls  the  matrix  inversion  subroutine  MINV. 

(g)  Computes  the  sum  of  the  squares  of  the  residuals 
between  the  calculated  and  experimental  data 
points . 

(h)   Cc?eputes  the  root-mean-square  residual  and  root- 
mean-square  error  (Reference  24)  of  the  current 
fit  to  date. 
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(i)   Computes  the  estimated  standard  deviations 

(Reference  24)  of  the  current  parameter  values. 

(j)   Tests  the  difference  between  the  current  and 
previous  values  of  the  root-mean-square  error 
to  determine  if  the  iteration  process  has 
converged. 

(k)   Computes  the  incremental  changes  for  the  para- 
meters (if  convergence  has  not  occurred). 

(1)   Computes  the  updated  parameter  values  (if 
convergence  has  not  occurred) . 

(m)   Returns  to  (c)  (if  convergence  has  not  occurred) . 

(n)   Computes  the  values  of  the  extracted  coefficients 
from  the  current  parameter  values  (after  con- 
vergence is  achieved) . 

(o)  Computes  the  estim.ated  standard  deviations  of 
the  extracted  aerodynamic  coefficients  (after 
convergence  is  achieved) . 

(p)   V/rites  extracted  aerodynamic  coefficients, 

estimated  standard  deviations,  and  the  pitch 
angle  output  that  represents  the  final  fit  to 
the  experimental  data  (after  convergence  is 
achieved) . 

Subroutine  functions 

The  names  and  functions  of  each  of  the  four  subroutines 
that  are  used  In  the  coefficient  extraction  program  are 
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given  below. 

ADDUM. --Subroutine  ADDUM  integrates  the  equation 
of  motion  and  parametric  differential  equations  using  a 
fourth  order  Runge-Kutta  method  for  starting  and  a  fourth 
order  Adams -Bashforth  predictor-corrector  method  for 
running.   It  calls  subroutines  XDOT  and  OUT.   This  sub- 
routine is  described  in  detail  in  Reference  26. 

XDOT . - -Subroutine  XDOT  computes  current  values  of 
first  derivatives  that  are  required  by  ADDUM. 

OUT. - -Subroutine  OUT  stores  the  results  of  the 
numerically  integrated  equation  of  motion  [Equation  (13)] 
and  parametric  differential  equations  [Equation  (15)]. 

MINV. - -Subroutine  MINV  inverts  an  K  x  K  matrix  using 
a  standard  Gauss- Jordan  technique  and  is  described  in 
detail  in  Reference  27. 

Required  input  data 

The  program  reads  six  categories  of  input  data. 
These  categories  and  the  specific  data  in  each  are  de- 
lineated in  Appendix  I .   The  format  and  units  of  the 
entries  on  a  specific  data  card  can  be  determined  from 
the  program  listing  and  nomenclature  list  provided  in 
this  appendix. 
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Program  optioas 

This  program  has  options  for  extracting  various 
combinations  of  aerodynamic  coefficients  from  the  given 
test  data,  in  addition  to  the  option  of  extracting  no 
coefficients  and  merely  integrating  the  equation  of  motion. 
These  options  are  controlled  by  the  numerical  value  of  the 
number  of  first  order  .equations  to  be  integrated,  N,  which 
is  read  by  the  program  on  the  first  data  card.   The  various 
options  are  giveB  in  Table  II. 

Table  II.   PROGRAiM  OPTIONS 

N  Coefficients  to  be  Extracted 


2   No  coefficients  are  extracted.   The  program 
integrates  the  equation  of  motion  for  the 
given  input  data  and  prints  the  results. 

10    "oj  '^Qy  ^mao'  ^mqo 

12    Og,  Qq,  Cjj^^q,  Cjj^qQ,  (^Yaa2 

14    Oq,  Oq,  C-^^q,    C^qQ,  C^^2>  ^mq2 

^^    °o'  "o>  ^mao'  ^mqo'  ^ma2'  '^mq2'  ^ma4 

^8    '^c*  "o'  ^maC  ^mqo'  '^ma2'  '^mq2'  ^ma4'  ^m5a 


Analysis  of  the  Sensitivity  to 
Noise  of  the  Chapman- Kirk  TechnTque 


The  evaluation  of  the  sensitivity  to  noise  of  the 
Chapman- Kirk  coefficient  extraction  technique  has  been 
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approached  from  several  directions.   Coefficients  have 
been  extracted  from  numerous  sets  of  computer  program 
generated  data  containing  only  measurement  errors  as  well 
as  data  containing  system  noise  and  measurement  errors. 

In  addition  to  the  use  of  artificially  generated  data, 
some  very  preliminary  work  has  been  done  with  noisy  wind 
tunnel  data  obtained  from  actual  experimentation.   The 
design  of  the  experiment  used  to  produce  the  dynamic 
data  is  described  by  Turner  in  Reference  28  along  with  some 
preliminary  results  obtained  with  the  one-degree-of -freedom 
program  described  earlier  in  this  chapter. 

Since  the  wind  tunnel  experimentation  is  still  in  a 
developmental  stage,  however,  the  results  presented  in  the 
remainder  of  this  chapter  are  those  obtained  from  arti- 
ficially generated  data. 

Generation  of  Noisy  Dynam.ic  Data 

The  computer  program  UFNOISE  described  in  detail  in 
Reference  29  was  used  to  generate  the  dynamic  data  from 
idiich  the  aerodynamic  coefficients  were  extracted.   This 
program  simulates  the  pitching  motion  of  symmetric  missile 
oscillating  in  a  wind  stream  about  a  fixed  point,  for  any 
given  initial  pitch  angle  displacement  and  initial  pitch 
angle  rate,  by  numerically  integrating  the  equation  of 
motion.   The  program  has  two  options  for  simulating  system 
noise:  it  considers  the  magnitude  and  direction  of  the 
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freestream  velocity  vector  as  normally  distributed  random 
variables,  with  programmer  set  mean  values  and  standard 
deviations,  to  determine  the  system  noise  perturbation 
accelerations  or  it  simply  selects  random  perturbation 
accelerations  which  have  zero  mean,  are  normally  distributed 
and  have  a  programmer  set  standard  deviation.   Both  methods 
are  essentially  equivalent.   New  values  of  the  system 
noise  perturbations  are  randomly  selected  at  each  numerical 
integration  step.   The  program  also  has  the  option  of 
simulating  random  measurement  errors  of  the  pitch  angle 
by  superimposing  iiormally  distributed  random  noise  on  the 
output  of  the  numerical  integration. 

The  basic  equation  of  motion  used  to  generate  the 
noisy  dynamic  data  was  a  slightly  simplified  form  of 
Equation  (13) : 

d  +  (C^  +  C6a2)a  +  (C3  +  C5a2)c(  =  wCt)  (16) 

subject  to 

a(0)  =  C,   ,     a(0)  =  C2    .  (17) 

The  true  values  of  the  physical  and  aerodynamic  con- 
stants used  in  generating  the  data  are  given  in  Table  HI. 
The  various  noise  level  standard  deviations  and  mean  values 
are  discussed  in  the  following  paragraphs  to  which  they 
are  pertinent. 
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Table  III.   TRUE  VALUES  OF  CONSTANTS 
USED  IN  GENERATING  DATA 


Constant  Value 

d(Ft)  0.333 

ACFt^)  0.0873 

ICsl'i.g  £t2)  0.1080 

V  (ft/sec)  500.0 

m 

q(lb/ft2}  297.0 

C^,Jrad-i)  -2.00 

C         Crad"3)  -24.5 

ma  2 

-60.0 
-163.0 
0.5235 
0.0 


c 

mqo 

%2f 

rad' 

^) 

a^Cra 

d) 

a^(ra 

d/sec) 

Effects  of  Random  Measurement  Errors 

Aerodyaaraic  coefficients  were  extracted  from  a  total 
of  nine  sets  of  data  containing  only  measurement  errors. 
Eacli  data  set  was  made  up  of  201  discrete  points  which  re- 
sult from  integrating  the  equation  of  motion  [Equation  (16)] 
in  increments;  of  0.005  second  for  a  total  of  1.000  second. 

The  desired  mean,  value  of  the  measurement  errors  for 
each  of  the  inine  data  sets  was  0.0  radian.   The  desired 
standard  devriations  of  the  random  errors  were  0.00146 
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radian,  0.00582  radian,  and  0,0174  5  radian,  with  three 
different  sets  of  data  being  generated  at  each  noise  level. 
These  standard  deviations  are  typical  of  measurement  un- 
certainty for  data  of  this  type  (Reference  9) . 

The  above  standard  deviations  correspond  to  percent 
noise  levels  of  1.1,  4.4,  and  13.2,  respectively  where 
percent  noise  level  is  defined  in  Reference  30  as 

■7 

Percent  Noise  = (18) 

Approximate  average  peak  amplitude 

with  a  being  the  standard  deviation  of  interest.   The 
approximate  average  peak  amplitude  can  be  determined  in  a 
variety  of  v.ays .   The  value  used  here  is  the  mean  positive 
peak  amplitude  for  the  three  cycles  that  result  when  the 
nominal  equation  of  motion  is  integrated  for  1.000  second. 

The  results  of  this  portion  of  the  analysis  are  given 
in  Table  IV  and  on  Figures  1  through  10. 

Table  IV  is  primarily  a  summary  of  the  percent  error 
in  the  extracted  aerodynamic  coefficients  together  with 
their  normalized  estimated  standard  deviations  for  the 
various  actual  measurement  errors,  a^.      The  percent  error 
in  the  extracted  coefficients  is  defined  by 


Percent  Error  (C.)  =  — ^ ^ (19) 

J  C. 
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where  C.  is  the  coefficient  estimate  determined  with  the 
program  in  the  fitting  of  the  data  and  C.  is  the  true 
value  of  the  coefficient  of  interest.   The  normalized 
estimated  standard  deviations  are  the  estimated  standard 
deviations  calculated  from  Equation  (12)  divided  by  the 
true  value  of  the  parameter  of  interest, 


Percent  a.      =  -J-  •  100      .  (20) 


Table  IV  also  contains  the  RMS  residual  for  each  fit  to  the 
noisy  data  and  the  percent  noise  levels  based  on  the  values 

An  exaiaple  of  one  of  the  program  fits  to  a  noisy  data 
set  is  shown  on  Figure  1.  Figures  2  through  5  show  the 
variation  of  the  percent  error  in  the  individual  extracted 
coefficients  with  measurement  error.  The  variation  of  RMS 
residual  with  measurement  error  is  shown  on  Figure  6.  The 
variations  of  the  normalized  estimated  standard  deviations 
of  the  extracted  coefficients  with  measurement  error  are 
shown  on  Figures  7  through  10. 

An  analysis  of  Table  IV  and  Figures  2  through  10 
reveals  the  following. 

(1)   Tiie  extracted  static  pitching  moment  coefficient 
derivatives  show  little  or  no  error  for  the 
entire  range  of  measurement  errors  considered. 
(See  Figures  2  and  3.) 
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(2)  Tfee  extracted  pitch  damping  coefficients  show 
some  error  for  the  lower  noise  values  (a  -0.005 
radian)  and  deviate  significantly  for  a  =0.018 
radian.   (See  Figures  4  and  5.) 

(3)  TIhe  variation  of  RMS  residual  with  measurement 
noise  is  essentially  linear  with  a  slope  of  unity 
(Figure  6).   This  indicates  that  the  RMS  residual, 
as  calculated  by  the  coefficient  extraction  pro- 
gra.m,  is  a  good  estimate  of  the  amount  of  measure- 
me^^at  noise  in  the  data  when  this  is  the  only 

type  of  noise  present. 

(4)  All  variations  of  the  normalized  estimated  standard 

deviations  of  the  extracted  coefficients,  a.     , 

'   jn' 

increase  linearly  with  measurement  error  (Figures 
7  through  10) .   Approximate  values  of  the  slopes 
o£  these  variations  are  given  in  Table  V.   These 
normalized  uncertainty  ratios  give  the  relative  (to 
eaich  other)  uncertainty  which  can  be  expected  when 
exitracting  coefficients  from  data  containing 
measurement  noise  using  the  given  coefficient 
extraction  program. 

(5)  T'.e  true  value  of  a  given  coefficient  is  contained 
within  the  interval  difined  by  its  extracted 
value  ±3c..  for  every  coefficient  extracted  from 
the  nine  sets  of  data. 
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Table  V.   EXTR.\CTED  COEFFICIENT  UNCERTAINTY  RATIOS 


Coefficient 


Approximate  Normalized 

Uncertainty  Ratio 

(Aa./Ao  ) 
*■   J    m^ 


C3 
C5 


2 

2 

5 

50 
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Effects  of  RaTidom  Measurement  Errors  and  System  Noise 

For  this  portion  of  the  analysis,  attempts  were  made 
to  extract  aerodynamic  coefficients  from  twelve  basic 
sets  of  dynainiic  data.   Each  data  set  was  made  up  of  201 
discrete  points  resulting  from  integrating  the  equation 
of  motion  in'  increments  of  0,005  second  for  4.00  seconds; 
the  pitch  angle  was  output  every  four  integration  steps 
so  that  the  itiine  between  data  points  was  0.020  second. 
The  time  betTs-seen  data  points  used  in  this  portion  of  the 
analysis  is  different  from  that  previously  used.   This 

change  was  msde  so  that  the  time  between  data  points  would 

2  8 
correspond  to  that  used  by  Turner    in  actual  experiment- 
ation. 

Of  the  twelve  basic  data  sets,  nine  contain  system 
noise  and  measurement  errors,  whereas  three  contain  only 
measurem.ent  errors.   Coefficients  v;ere  extracted  from  these 
latter  three  to  determine  if  the  above-mentioned  change  in 
the  time  increment  between  data  points  had  any  appreciable 
effect  on  the  extracted  coefficients  or  their  uncertainties 
The  estimates  of  the  uncertainties  were  improved  due  pri- 
marily to  the  fact  that  a  longer  data  record  was  being 
analyzed. 

The  desired  mean  values  of  the  freestream  velocity 
vector  direction  fluctuations  and  measurement  errors  for 
all  data  sets  were  zero.   The  desired  standard  deviation 
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of  the  measurement  errors,  o   was  0.00582  radian 
(3o   --  1.0  degree)  for  all  sets;  the  desired  standard 
deviation  of  the  velocity  magnitude  fluctuations,  a 
was  5,0  ft/sec  for  all  nine  data  sets  containing  system 
noise.   The  desired  standard  deviations  of  the  velocity 
direction,  c-  ,  were  0.02910  radian  (3a   =  5  degrees), 
0.05820  radian  (3o^  =  10  degrees),  and  0.11640  radian 
(3o   =  20  degrees) ,  with  three  different  sets  of  data 
being  generated  at  each  noise  level. 

The  random  velocity  fluctuations  act  as  a  forcing 
function  for  the  equation  of  motion  and  cause  oscillations 
even  if  the  vehicle  has  no  initial  displacement.   The 
resultant  maximum  amplitudes  of  these  forced  oscillations 
are  usually  within  a  certain  magnitude  or  noise  band  width. 
The  widths  of  the  noise  bands  for  the  data  used  in  this 
analysis  are  generally  equivalent  to  the  3a   values,  and 
these  have  been  used  to  calculate  the  percent  system  noise 
values.   The  approximate  average  peak  amplitude  used  in  tlie 
percent  noise  calculations  is  the  mean  positive  peak  ampli- 
tude for  the  10  cycles  that  result  when  the  nominal  equation 
of  motion  is  integrated  for  approximately  four  seconds. 

The  results  of  tliis  portion  of  the  analysis  are  given 
in  Table  VI  and  on  Figures  11  through  19. 

Table  VI  is  a  summary  of  the  percent  error  in  the 
extracted  aerodynamic  coefficients  and  their  normalized 
estimated  standard  deviations  along  with  the  various  noise 
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Table  VI.   Smi>IARY  OF  RESULTS  FOR  EXTRACTING  AERODYNAMIC 
(COEFFICIENTS  FROM  DYNAMIC  DATA  CONTAINING 
RANDOM  MEASUREMENT  ERRORS  AND  SYSTEM  NOISE 


Data 
Set 

(Had) 

(Ft/sec) 

^a 
(Rad) 

Percent 

System 

Noise 

10 

0.O0527 

0.00 

0.00 

0.00 

11 

0.00589 

0.00 

0.00 

0.00 

12 

0.OT574 

0.00 

0.00 

0.00 

13 

0.00611 

5.05 

0.02i     J 

41.4 

14 

0.00597 

5.18 

0.03038 

4  3.1 

15 

0,00587 

4.88 

0.02828 

40.2 

16 

0.OO603 

5.01 

0.005516 

78.4 

17 

0.00  562 

4.90 

0.06120 

8  6.8 

18 

0.O0600 

4.96 

0.05662 

SO. 3 

19 

0.O0S72 

5.08 

0.1126  2 

152.9 

20 

0.00547 

4.85 

0.11706 

1S8.9 

21 

0„O0571 

4.96 

0.12058 

163.7 
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Table  VI.   (Extended) 


Percent  Error  in  Extracted  Coefficients  and  RMS 

Normalized  Estimated  Standard  Deviations         Residual 


03(03)         Ci+(a4)  05(05)  C6(a6) 


(Rad) 


-0.27(0.25)   -0.84(0.93)    0.68(0.65)  5.72(12.92)  0.00520 

-0.34(0.28)    0.29(1.04)    1.51(0.73)  -3.60(14,43)  0.00574 

0.42(0.27)   -0.36(1.01)   -0.94(0.69)  4.40(13.88)  0.00563 

0.55(0.64)    8.25(2.75)    9.55(1.88)  329.59(43.58)  0.01253 

4.19(0.63)    5.84(2.41)    1.35(1.63)  -44.80(34.40)  0.01244 

2.87(0.55)   20-73(2.10)   -2.80(1.24)  -247.20(26.39)  0.01026 

2.71(1.24)  -42.80(3.95)    6.54(3.64)  605.77(66.79)  0.02806 

14.44(1.32)   -0,90(5.35)    5.44(3.87)  197.10(90.61)  0.02627 

1.41(0.70)  -35.67(2.33)   14.59(2.30)  588.34(44.19)  0.01488 

DIVERGED 
DIVERGED 
DIVERGED 
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level  standsxd  deviations,  percent  system  noise  levels, 
and  the  RiMS  residual  of  each  fit  to  a  noisy  data  set.   The 
percent  system  noise  levels  are  similar  to  those  encountered 
experimentally  by  Turner  in  Reference  28. 

Figure  11  is  an  example  of  one  of  the  program  fits 
to  a  noisy  data  set.   Figures  12  through  15  show  the  varia- 
tion of  perG:ent  error  for  the  extracted  coefficients  with 
system  noise.   The  variations  of  the  normalized  estimated 
standard  deviations  of  the  extracted  coefficients  with 
system  noise  level  are  shown  on  Figures  16  through  19. 

An  analysis  of  Table  VI  and  Figures  12  through  19 
reveals  the  following: 

(1)  Tbe  errors  in  the  linear  sti'tic  pitching  moment 
coefficient  derivative  parameter,  C3,  are  gen- 
erally less  than  5  percent  for  the  entire  range 
o£  system  noise  considered.   The  errors  in  the 
corresponding  nonlinear  term,  C5,  are  generally 
less  than  10  percent.   (See  Figures  12  and  13.) 

(2)  TIrs'e  pitch  damping  parameters  show  significant 
error  for  all  nonzero  system  noise  levels.   The 
extracted  values  of  the  nonlinear  term,  Cg ,  are 
sometimes  in  error  by  several  orders  of  magnitude 
at  the  higher  noise  levels.   (See  Figures  14 

and   15.) 

(3)  'The   estimated  standard  deviations  of  the  extracted 
coefficients  generally  increase  as  the  system  noise 
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increases.   The  estimated  standard  deviations 
of  coefficients  extracted  from  dynamic  data 
containing  both  system  noise  and  measurement 
errors  are  generally  too  small  and  do  not  reflect 
the  true  uncertainty  of  the  extracted  coeffi- 
cients as  was  true  in  the  previous  case  when 
tSue  data  contained  only  measurement  errors. 
(See  Table  VI  and  Figures  16  through  19.) 

(4)  All  attempts  at  extracting  the  complete  set  of 
coefficients  from  data  with  a  system  noise  band 
of  approximately  20  degrees  failed  (see  Table  VI) , 
The  failures  resulted  when  a  value  of  Cg  was 
e\'-entually  calculated  by  tho  iterative  process 
wkich  caused  the  solution  to  the  equation  of 
notion  to  diverge. 

(5)  The  divergence  problem  can  sometimes  be  circum- 
vented by  not  attempting  to  extract  C    from 
extremely  noisy  data.   The  resulting  coefficients 

that  are  extracted,  C    ,  C   „,  and  C    ,  have 

m.ao    ma  2'       mqo 

accuracies  comparable  to  those  extracted  from 
noisy  data  with  a  ±10  degrees  system  noise  band. 
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CHAPTER  IV 

DEVELOPMENT  OF  THE  EXTENDED  KALMAN  FILTER  FOR 
ESTIMATING  PARAMETERS  AND  THEIR  UNCERTAINTIES 


This  chapter  presents  an  abridged  derivation  of  the 
Kalman  filter  for  discrete  and  continuous  linear  systems, 
followed  by  a  general  statement  of  the  extended  Kalman 
filter  for  continuous  nonlinear  systems.   The  extended 
filter  is  then  used  as  a  base  for  the  development  of  a 
specific  algorithm  for  estimating  states  and  parameters 
of  the  second  order  equation  of  pitching  motion  for  a 
missile  being  forced  by  random  disturbances. 

Development  of  the  Kalman  Filter 

The  original  derivation  of  the  Kalman  filter  was 

I  Q 

presented  by  Kalman    in  1960  for  multistage  systems 
making  discrete  linear  transitions  from  one  stage  to 
another.   Kalman  and  Bucy   gave  the  analogous  development 
for  continuous  linear  ..ystems  approximately  one  year  after 
the  first  work  was  published.   The  purpose  of  the  linear 
filter  is  to  provide  estimates  of  the  state  of  a  system 
hy   making  use  of  measurements  of  all,  or  some,  of  the 
state  vector  components  of  the  system.   The  system  is 
assumed  to  be  operating  in  the  presence  of  random  distur- 
bances, the  statistical  properties  (i.e.,  mean  and  variance) 
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o£  which  are  known.   The  measurements  of  the  state  vector 
components  of  the  system  are  also  assumed  to  have  random 
errors  of  known  statistics. 

In  addition  to  the  original  derivations  of  Kalman 
and  Bucy,  a  number  of  other  methods  offering  various  degrees 

of  insight  but  leading  to  the  same  results  are  available. 

2  0 
A  brief  development  taken  primarily  from  Bryson  and  Ho 

is  presented  here.   Other  derivations  or  developments  of 

31         32 
the  filter  equations  are  given  by  Jazwinsky,    Sorenson 

and  Barham. 

The  treatment  given  here  starts  with  a  static  system, 

is  extended  for  a  single-stage  linear  transition  which 

leads  directly  to  linear  multistage  processes.   Finally, 

by  making  use  of  a  limiting  process,  the  desired  form  of 

the  equations  for  a  continuous  linear  dynamic  system  is 

given. 

Static  Systejn 

The  problem  at  hand  is  to  estimate  the  n-component 
state  vector  X  of  a  static  system  using  the  p-component 
measurement  vector,  z,  containing  random  errors,  v,  which 
are  independent  of  the  state.   The  measurement  vector  can 
be  represented  as 

z  =  KX  +  V       ,  (21) 
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where  H  is  a  known  p  x  n  matrix.   Conditions  on  the  measure- 
ment errors  are 

E(v)  =  0  (22) 

E(vv'^)  =  R     ,  (23) 

where  R  is  a  knovm  matrix  of  dimension  p.   It  is  assumed 
that  a  prior  estimate  of  the  state,  designated  as  X,  is 
available  and  also  that  the  covariance  of  the  prior  esti- 
mate, M,  is  known.   Thus 

M  =  E[(X-X)(X-X)'^]      ,  (24) 

where  M  is  of  dimension  n. 

The  desired  estimate  of  X,  taking  into  account  the 
measurement,  z,  is  the  weighted- least-squares  estimate, 
X,  which  minimizes 

J  =  [(X-X)"^  M-i  (XX)  +  (z-HX)"^  R-i  (z-lIX)].       (25) 

Differentiating  Equation  (25)  with  respect  to  X,  setting 
the  resulting  equation  to  zero  and  solving  for  X  yields 

X  =  X  ••■  Ph'^R-1  (z-HX)      ,  (26) 

where 

p-i  =  M-l  H  H-^R-Hl         .  (27) 
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The  quantity  P  is  the  covariance  matrix  of  the  error  in 
the  state  estimate  after  measurement,  X,  i.e., 

P  =  E[(X-X3(X-X)'^]      .  (28) 

Single -Stage  Transitions 

It  is  now  desired  to  estimate  the  state  of  a  system 
that  makes  a  discrete  transition  from  state  0  to  state  1 
according  to  the  linear  relation 

^1  --    ^^o  "  ^0%     '  .C29) 

where  <})   is  a  known  transition  matrix  of  dimension  n  and 
r   is  a  knovrn  n  x  r  matrix.   The  mean  and  variance  of  the 
random  forcing  vector  are  assumed  to  be  known  and  are 
given  by 

E(w^)  =  w„    ,     E[Cw^-w  )  (w  -w  )'^]  =  Q     .   (30) 

^O"^      O     *        "-^OOOO         ^O         V"-/ 

The  statistical  properties  of  the  random  state  vector  are 
assumed  to  be  known  initially  as 

ECX^)  =  i^    ,     E[(X^-X^)(i^-X^)T]  =  p^    .   (31) 

It  is  also  assumed  that  X^  and  w  are  independent.   From 

o  o  ^ 

this    information  Xi    is    a   random  vector   vvhose   mean  value 
and   covariance    are 

X,    =    $^i      +    r   w^  (32) 

loooo  ^      ^ 
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M,  =  (j)  P  4>'^  +  r^Q  r^     .  (33) 

1     o  o  o     o^o  o  ^      ■' 


Suppose  now  that  measurements  of  the  state  are  made  after 
the  transition  to  state  1;  then  from  Equations  (26)  and 
(27)  the  best  estimate  of  X   is 

Xj  =  Xj  +  PiH'[r-i  (z-H^X^)     ,  (34) 

where 

Pj  =  (M'^  +  h'[r-1h^)-i 

=  Ml  -  MjH'J  (H^M^h]"  +  Rj^)  HjMi      .       (35) 

Linear  Multistage  Processes 

The  developments  of  the  two  preceding  sections  can 
be  extended  for  linear,  stochasitc,  multistage  processes. 
Given  the  following  difference  equation  system  model 

X.         =    (J).X.    +    r.w.  ,  (36) 

i+iiiii  *  ^      ■^ 

where 


E(X^)    =   X^  (37) 

ECw.)    =   w^  (38) 

E[(X    -X    )  (X    -X^)*^]  =   M^                                                          (39) 

'^^OOOO'^  o 
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E[(w^-w.)(w^-w^.)]  =  Qi«ij  (40) 

E[(w.-w.)CXq-X^)'^]  =  0  (41) 

and  measurements  o£  the  state 

z.  =  H-X.  +  V.        ,  (42) 

where 

'< 

E(v^)  =  0  (43) 

E(v.vJ)  =  R.6.^  (44) 

E[(w.-M^)v]']  =  0     ,     E[(X^-X^)v^]  =  0    ,      (45) 

the  estimate  of  the  state  is 

X.    =   X.    -H   K-(z.-H.X.)                ,  (46) 

where 

.^i.i    =    *i^i    ^    ^i^i  (47) 

Ki    =    PiHiRi"'  (43) 

P.    =     (M.-l    +    hTr. -iH.)'^ 
1^1  111-' 

-    M.     -    M.H't'(H.M.HT    +    R.-1)H.M.  (49) 

1  11^111  I'^ii  ^      ■' 


^'i.i  =  h^A  '  'i^A       '  (50) 
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Equations  (46),  (47)  and  (48)  are  the  discrete  Kalman 
filter  with  the  state  variable  covariances  given  by- 
Equations  (49)  and  (50) .   As  can  be  seen  from  the  above 
equations,  the  Kalman  filter  is  essentially  the  same  as 
the  system  model  [Equation  (36)].   The  differences  are 
(1)  the  actual  system  noise,  which  is  random  from  one  stage 
to  the  next  in  Equation  (36)  ,  is  replaced  by  its  mean  or 
expected  value,  and  (2)  there  is  a  correction  term  based 
on  the  difference  between  the  actual  measurement  of  the 
state  and  its  predicted  value.   The  difference  term  is 
multiplied  by  a  gain,  K^,  which  is  essentially  the  ratio 
of  the  uncertainty  in  the  state  to  the  uncertainty  in 
the  measuremeat.   If  the  covariances  of  the  measurement 
errors  are  large,  the  gain  will  be  relatively  small  and 
the  corresponding  difference  term  will  have  little  effect 
on  the  estimate  of  the  state;  if,  however,  the  system 
noise  is  relatively  large  or  the  measurement  errors  are 
very  small,  the"  gain  will  be  large  and  differences  between 
the  actual  and  predicted  measurements  of  the  state  at  a 
given  stage  have  increased  significance  in  the  state 
variable  estimates. 

The  Continuous  Kalman  Filter 

By  applying  a  limiting  process  to  the  discrete  filter 
with  the  tine  between  stages  tending  to  zero,  the  linear 
system  model  becomes 


X  =  F(t)X  +  G(t)w(t) 


and  the  continuous  Kalman  filter  is  given  by 


X  =  FX  +  Gw  +  PH^R-i (z-HX) 

f  =  FP  +  PF^  +  GQg'^  -  Ph'^R-IHP 


(51) 

(52) 
(53) 


The  Extended  Kalman  Filter 


The  extension  o£  the  Kalman  filter  for  estimating 
the  states  of  nonlinear  systems  in  the  presence  of  noise 
has  been  given  by  Bryson  and  Ho    and  Jazwinski,"   among 


others,  as 

X  =  f(X,t)  +  G(t)w(t)  +  P 


I 

II]   R-i[z(t)-h(X,t)]    (54) 


3f 


P  +  P 


if] 


GQG'  -  P 


'9h]'  R-l|lil" 

axj  ^  [ax. 


for  the  nonlinear  system 


X  -  f(X,t)  +  G(t)w(t) 


with  measurements 


wnere 


z(t}  -  h(X,t)  +  v(t) 


F,[wrt)l  =  w(t) 


■e[w(t)-w(t)][wCf)]  =  Q(t)d(t-t') 


E[X(t  Jj  -  X 
*-  ^  o  ■■    o 


E[x(t^)-x  j[xrt,)-x^j'  -  p, 

'•  ^  o'   0-'  ■    (j-^   o-"      o 


(55) 


(56) 

(57) 

(58) 
(59) 
(60) 
(61) 
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E[X(tQ)-X^][w(t)-w(t)]'^  =  0  (62) 

E[v(t)]  =  0  (63) 

E[v(t)v'^(t' )]  =  R(t)6(t-t')  (64) 

E[X(t^)-X^][v(t)]'^  =  0  (65) 

E[w(t)-vv(t)]  [v(t' )]'^  =  0  .            (66) 


The  Parameter  Estimation  Algorithm 

In  this  section  a  specific  parameter  and  state  esti- 
mation algorithm  using  the  extended  Kalman  filter  is 
developed  for  the  basic  equation  of  pitching  motion  pre- 
viously analyzed  with  the  Chapman-Kirk  technique. 

The  nonlinear  equation  of  motion,  or  system  model, 
is  tlius 

a    +     (C^+Cga^)^  +  (C3+C5a2)a  =  w(t)  (67) 

where 

E[w(t)]  =  0    ,     E[w(t)w(t')]  =  q5(t-t')     .    (68) 

The  measure?r,8nt5  of  the  state  of  the  system  are  assumed 
to  be  given  by 

zrt)  =  a(t)  +  v(t)  (69) 


waere 


Etv(t)]  -  0     ,     E[v(t)v'^(t')]  =  r5(t-t')    .   (70) 
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Now  to  reduce  Equation  (67)  to  the  required  system  of 
first  order  differential  equations,  let 


a  =  X2         .  (71) 

At  this  point,  to  estimate  the  aerodynamic  parameters 
appearing  in  the  equation  of  motion  in  addition  to  the 
state  variables,  the  state  vector  is  augmented  by  setting 

C5  =  X5 

Ce  =  Xg  (72) 

with  the  constraints 

X3  =  0 

Xu  =  0 

X5  -  0 

Xs  =  0         .  (73) 

Making  use  of  Equations  (71) ,  (72)  and  (73)  ,  Equation  (67) 
now  reduces  to  the  following  nonlinear  system  of  first 
order  equations: 
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Xi    = 

X 

X2     = 

- 

X3  = 

0 

X.    = 

0 

Xs    = 

0 

X5    = 

0 

(X,+X5Xi2)X2-CX3+X5Xi2)Xi  +  W(t) 


(74) 


with  linear  laeasurements 


=  Xi  +  V 


(75) 


Comparing  Equations  (74)  and  (75)  with  (56)  and  (57)  the 
following  mattrices  may  be  identified: 


G(t) 


H(t)  = 


f(X,t) 


X2 

-(X4+XsXi2)X2 

0 
0 
0 
0 


-   (X3+X5Xi2)Xi 


(76) 


Now,  since  the  mean  value  of  the  system  noise  is  assumed 
to  be  zero  ^nquation  (68)]  and  because  the  measurements 
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of  the  state  are  related  linearly  to  the  state  [Equation 
(75)],  the  extended  Kalman  filter,  previously  given  by 
Equations  (.54^    and  (55)  ,  can  be  simplified  somewhat  to 


X  =    (X,t)    +   PH^    -    [z-HX] 


P   =  —  P   + 
^         8X   ^ 


ph'^ 

i' 

r 

V        J 

p 

hf] 

,9Xj 

T 

T  T 

+    GqG       -    PH^ 


.r. 


HP 


(77) 
(78) 


wliere 


H=[l      0      0      0      0      0] 


Applying  Equations  (76)  and  (77)  yields 


• 

Xi 

• 

L 

• 

X3 

• 

X, 

X5 

^6 

>.        ; 

(X3+X5Xi2)Xi 


X2 
(X4+Ar;Xi-)X2 

0 
0 
0 


Pii  P12  Pi3  Pi^  Pi5  Pig 


Pi  2  ?2  2  P23  P 

Pl3  P23  P33  P 


24  -^25  ^26 


31+  Ms 


P^, 


Pl4     ^^2h     P31+     Pi+it     Pi+5     Pit6 


Pl5  P25  P35  P45  P55  P56 

Pi6  ?2S  Psb  P^e  P56  Pes 


1 

0 
0 
0 
0 
0 

I  J 


:(z,-  [10  0  0  0  0] 


Xi 

X2 

X3 

X4 

X5 

Xe 
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Carrying  out  the  prescribed  multiplications  and  equating 
like  components  of  the  resulting  two  column  matrices 
yields  the  desired  filter  equations, 

i     .     Pi  1 

X2  =  -(i^+XgX2)X2  -  (X3+X5x2)Xi  +  -^   Cz-Xi) 

•  Pl3 

X,  =  ^r'-    (z-X,) 

•  Ft  s     - 
X5  =  -^  (z-X,) 

Xg  =  -l-^   (z-Xi)  .  (79) 


The  necessary  covariance  equations  are  obtained  from  the 
matrix  Ricatti  equation  [Equation  (78)]  making  use  of 
Equations  (76)  and  the  fact  that 


o  r 


3X 


X=X 


- 2X1X2X5-X3-3X1X5 
0 
0 
0 
0 


1 

0 

0 

0 

0 

(X,+X6 

-   9 

X'l) 

-Xi 

-X2 

-Xi 

-X1X2 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

(80) 
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The  results  o£  substituting  Equations  (76)  and  (80)  into 
(78) ,  carrying  out  the  rather  tedious  matrix  multipli- 
cations and  equating  like  terms  are 


Pii 

Pi  2 
Pl3 
Pl^ 
Pl5 
Pl6 

• 

P22 
P23 
P24 
P25 
P26 
P3  3 
P3.. 
P35 
P36 
P^U 
P45 

P^e 
P55 
P56 
Pee 


2Pi2-^Pii' 

P22-APii-BPi2-XiPi3-X2Pm-CPi5-DPie4PiiPi2 

P234P11P13 

P2t,4^iiPi't 
P25-7P11P15 
P26-7PiiPie 

q-2(APx2  +  BP22  +  XiP2  3  +  X2P24-^CP2  5  +  nP2  6)-7Pl2^ 
-APi3-BP2  3-XiP3  3-X2P34-CP3  5-DP3e-^Pi2Pl3 


-AP 

-AP 
-AP 

It, 


■Ip 
r 


-ip 
r 


-¥ 


r 

h 
r 

1 


r 

Ir 


■BP 


2'+"^lP3i+'X2P4'+ 


■CP 


4b' 


DP.A--P 


«f6 


12^14 


5-BP2  5-XiP3  5-X2P4  5-CP55-0P56-fPl2Pl5 


6-BP26 


3 
3P14 


3P1  5 


3^16 

2 

k 

.'.Pi  5 

^Pl6 

2 

5 

5P15 
2 


XiP,.-XoP, 


l''3  6-^2^'+6-CP5  6-^^66-T^12^16 


DP. 


(81) 


where 
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A  =  2X1X2X5  +  X3  +  3X1X5 


B  =  X^  >  XgXf 

C  =  Xi 

B  =  X1X2  .  (82) 

The  desired  estimates  of  the  states  and  parameters 
are  obtained  by  numerically  integrating  the  filter  equations 
[Equations  (79)]  and  their  covariances  [Equations  (81)]. 
Initial  estimates  of  the  states  and  their  covariances  as 
v\fell  as  the  variances  of  the  system  noise  and  measurement 
errors  are  assumed  to  be  available.  In   the  event  that  the 
data  consist  of  discrete  measurements,  which  is  usually 
the  case,  the  constant  time  between  data  points  should  be 
equal  to  the  numerical  integration  step  size.   This  is 
necessary  because  the  state  estimates  are  updated  at  each 
integration  increment  and  in  doing  this  the  filter  requires 
knowledge  of  the  difference  between  measured  and  estimated 
state  values. 

If  the  filter  successfully  adjusts  the  true  states  and 
the  imbedded  parameters  so  that  Xi  is  an  adequate  match  to 
the  data,  z,    the  time  derivatives  of  the  filter  equations 
for  t!ie  pararaeters  become  small  [see  Equations  (79)]  and 
the  parameters  reacli  steady-state  or  near  steady-state 
conditions . 


CHAPTER  V 

USE  OF  THE  EXTENDED  KALMAN  FILTER  FOR 
ESTIMATING  PARAMETERS  AND  THEIR  UNCERTAINTIES 


This  chapter  presents  an  analysis  using  the  extended 
Kalman  filter  with  parameter  augmented  state  vector  to 
determine  aerodynamic  coefficients  and  their  uncertainties 
from  noisy  dynamic  data.   A  description  of  the  digital 
computer  program  used  in  the  analysis  is  given  first, 
followed  by  z^esults  obtained  with  this  program  from  a 
variety  of  noisy  data  sets. 

The  Extended  Kalman  Filter  Program 

The  extended  Kalman  filter  program  is  written  in 
Fortran  TV  for  use  primarily  on  an  IBM  model  360/65 
digital  computer.   The  basic  function  of  the  program  is 
to  integrate  numerically  the  27  first  order  differential 
equations  which  are  given  in  the  previous  chapter,  and 
which  comprise  the  extended  Kalman  filter  with  parameter 
augmented  state  vector.   The  specific  functions  of  the 
main  program  and  its  subroutines,  as  well  as  the  required 
input  data,  are  described  in  the  following  paragraphs. 
A  listing  of  the  complete  program  is  given  in  Appendix  II. 
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Main  Program  Functions 

The  main  program  reads  and  writes  the  input  data 
and  calls  the  numerical  integration  subroutine,  ADDUM. 

Subroutine  Functions 

The  names  and  functions  of  each  of  the  three  sub- 
routines that  are  used  in  the  extended  Kalman  filter 
program  are  given  below. 
ADDUM 

Subroutine  ADDUM  integrates  the  27  filter  and 
covariance  equations  using  a  fourth  order  Runge-Kutta 
method  for  starting,  followed  by  a  fourth  order  Adams- 
Bashforth  predictor-corrector  method  for  running.   It 
is  essentially  identical  to  the  subroutine  of  the  same 
name  used  in  the  Chapman-Kirk  coefficient  extraction 
program  and  is  described  in  detail  in  Reference  26. 
Subroutines  XDOT  and  OUT  are  called  from  ADDUM. 
XDOT  .       . 

Subroucine  XDOT  computes  values  of  tiie  time  deri- 
vatives for  the  differential  equations  being  integrated 
by  ADDUM. 
OUT 

This  subroutine  writes  the  output  of  the  numerical 
integration. 
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Required  Input  Data 

The  program  reads  four  categories  o£  input  data.   These 
categories  and  the  specific  data  in  each  are  delineated  in 
Appendix  II.   The  format  and  units  of  the  entries  on  a 
specific  data  card  can  be  determined  from  the  program 
listing  and  nomenclature  list,  which  are  also  given  in 
Appendix  II. 

Analysis  of  the  Extended  Kalman  Filter 


The  use  of  the  extended  Kalman  filter  for  estimating 
the  parameters  of  interest  is  analyzed  for  linear  and  non- 
linear systems  v;itli  data  containing  only  measurement  errors 
as  well  as  data  containing  botli  system  noise  and  measurement 
errors.   The  analysis  begins  with  a  linear  system  and  data 
containing  only  measurement  errors  and  progresses  through 
increasing  stages  of  difficulty  up  to  nonlinear  systems  and 
data  v\/hich  have  been  corrupted  by  both  system  noise  and 
measurement  errors.   The  previously  discussed  computer 
program  and  algorithm  require  no  modifications  to  consider 
the  linear  case;  by  initially  setting  the  parameters  that 
are  the  numerical  coefficients  of  the  nonlinear  terms  in 
the  equation  of  motion  (and  their  variances)  equal  to  zero, 
the  extended  filter  for  a  nonlinear  system  reduces  to  the 
one  required  for  a  linear  system. 
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All  data  used  in  the  analysis  were  generated  by  the 
previously  mentioned  computer  program,  UFNOISE  (Reference 
29) .   The  true  values  o£  the  constants  used  in  the  data 
generation  are  those  previously  given  in  Table  III  with  the 
exception  of  the  initial  value  of  the  pitch  angle  for  the 
linear  system.   This  initial  displacement  is  a  more  realistic 
0.1745  radian  for  the  linear  cases  considered. 

Linear  Systems  with  Measurement  Errors  Only 

After  the  extended  Kalman  filter  computer  program  had 
been  constructed  and  checked,  it  became  apparent  that  the 
technique  would  probably  be  best  understood  by  considering 
relatively  simple  cases  at  first  and  then  progressing  to 
more  difficult  ones  as  confidence  in  the  technique  was  gained. 
To  this  end,  the  first  success  with  the  program  was  realized 
when  analyzing  data  containing  only  measurement  errors  for 
a  linear  system. 

The  pitch  angle  data  used  in  this  part  of  the  analysis 
were  generated  by  integrating  the  equation  of  motion 

a  +  C^i  +  Cja  =  0  (83) 

and  superimposing  random  Gaussian  errors  on  the  output. 
The  standard  deviation  of  the  errors  is  0.00588  radian, 
which  corresponds  to  3o  measurement  errors  of  approximately 
1.0  degree.   The  numerical  integration  step  size  used  in 
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generating  the  data  was  0.005  second  and  the  pitch  angle 
was  output  every  integration  increment.   The  equation  was 
integrated  for  a  total  of  2.00  seconds. 
Basic  filter  performance 

The  results  of  using  the  extended  Kalman  filter  to 
identify  the  correct  values  of  the  parameters  of  interest, 
X3  and  Xi+,  are  shown  on  Figures  20  and  21.   These  figures 
show  time  histories  of  the  percent  errors  in  the  estimated 
values  of  the  parameters  that  are  computed  by  the  extended 
filter  program  when  analyzing  the  noisy  data  described 
previously.   The  percent  error  is  defined  by 

(X.-X.).IOO 

PERCENT  ERROR  (X.)  =  — ^ ^ '  (84) 

^         X. 

The  input  for  the  measurement  error  variance  was 
the  one  previously  computed  by  the  simulation  program 
based  on  the  errors  that  were  actually  put  into  the  data. 
The  initial  values  of  the  parameter  variances  were  computed 
from  knowledge  of  the  true  value  of  the  parameters  and  the 
arbitrarily  selected  initial  values  of  the  parameters. 
These  sample  variances  are  defined  by 

P.,(0)  =    [X.(0)-X  (0)]2  .  (85) 

J  J  J  J 
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This  method  of  initializing  the  pai'ameter  variances  was 
chosen  merely  to  generate  the  required  initial  numerical 
values.   In  actual  experimentation,  the  initial  variance 
values  would  depend  on  the  method  of  selecting  the  initial 
parameter  estimates  and  prior  knowledge  as  to  how  accurately 
these  initial  parameter  estimates  were  in  relation  to  the 
true  values  of  the  parameters. 

The  initial  values  of  all  covariances  were  chosen  as 
zero.   This  implies  that  errors  in  the  estimates  of  the 
individual  state  vector  components  are  initially  uncor- 
related.   The  initial  values  of  the  variances  for  the 
pitch  angle  and  pitch  angle  rate  were  also  chosen  as  zero. 
For  wind  tunnel  dynamic  experimentation  where  the  model 
is  initially  held  rigid  at  some  given  displacement  to  the 
wind  stream,  this  seems  to  be  a  valid  assumption. 

As  can  be  seen  from  Figures  20  and  21,  the  filter 

does  an  excellent  job  of  identifying  the  parameters  of 

interest.   The  X3  parameter,  wliich  corresponds  to  the 

C   term,  and  which  is  initially  in  error  by  25  percent, 
ma      '  ^  J  V  i 

is  identified  with  approximately  zero  percent  error  in 
less  than  0.3  second.   The  correct  identification  of  the 
damping  paiameter  takes  slightly  more  time  but  the  results 
are  of  equal  quality. 
Effects  of  variations  in  the  initial  parameter  variances 

This  section  presents  some  effects  on  the  near  steady- 
state  estimates  of  the  parameters  of  interest  and  their 
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near  steady-state  variances  for  various  values  o£  the 
initial  parameter  variances.   The  parameter  variances 
were  initialized  to  values  that  were  25  percent  higher 
than  those  used  in  the  previous  section  and  also  to 
values  that  were  25  percent  lower  than  the  referenced 
values.   The  effects  of  these  initializations  are  shown 
on  Figures  22  through  25  and  in  Table  VII  for  the  linear 
system  of  interest. 

Table  VII  is  a  summary  of  the  percent  errors  in  the 
parameter  variances  and  the  normalized  parameter  uncertain- 
ties for  the  above  mentioned  variations  in  the  initial 
parameter  variances. 

The  same  results  are  depicted  graphically  on  Figures 
22  and  23.   The  numerical  results  in  Table  VII  are  based 
on  the  near  steady-state  parameter  and  variance  values 
that  result  after  the  filter  has  integrated  for  1.5 
seconds.   These  correspond  to  the  last  point  shown  on  the 
supporting  figures.   The  normalized  parameter  uncertainties 
are  defined  by 


/P77 
PERCENT  (/PT~)   =  — li  •  100      .  (86) 

JJ  n     ^j 


Time  histories  of  the  parameter  variances,  P33  and  Pt+u, 
are  shov\rn  on  Figures  24  and  25.  All  of  the  information 
presented  on  these  figures  indicates  that  initial  values 
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Table  VII.   EFFECTS  OF  VARIATIONS  IN  THE  INITIAL 
PARAMETER  VARIANCES  ON  NEAR  STEADY- 
STATE  PARAMETERS  AND  THEIR  VARIANCES 
(LINEAR  SYSTEM,  MEASUREMENT  ERRORS  ONLY) 


Variation  in      Errors  in  Near 
Initial  Parameter     Steady-State 


Normalized  Near 

Steady-State 

Variances 

(Percent) 

v^P33        ^'P^k 

n         n 

Variance*  Parameters 

(Percent)  (Percent) 

X3  X^^ 


25  0.17        0.55        0.12       0.91 

0  0.19        0.46        0.11       0.91 

25  0.22        0.38        0.11       0.91 


relative  to  sample  variance 
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o£  parameter  variances  in  the  range  investigated  have  no 
effect  on  final  parameter  estimates  or  their  variances 
when  using  the  extended  Kalman  filter  to  analyze  data  for 
a  linear  system  containing  only  measurement  errors. 

Linear  Systems  with  Measurement  Errors  and  System  Noise 

Data  for  this  portion  of  the  analysis  were  generated 
with  the  computer  program  UFNOISE  by  numerically  integrating 
the  following  equation  of  motion, 

a  +  C^i  +  Cga  =  w(t)  (87) 

and  superimposing  random  Gaussian  errors  on  the  output. 
The  forcing  function,  w(t),  is  also  random  in  nature  and 
Gaussian  wi£h  zero  mean.   The  standard  deviation  of  the 
system  noise  was  4.80  rad/sec^,  which  is  of  sufficient 
magnitude  to  drive  the  oscillations  in  a  noise  band  of 
approximately  0.0872  radian  (5  degrees)  for  zero  initial 
displacement.   The  standard  deviation  of  the  measurement 
errors  was  0.00556  radian.   Also,  as  before,  the  equation 
of  motion  was  integrated  in  increments  of  0.005  second  for 
a  total  of  2.00  seconds. 
Basic  filter  performance 

The  basic  performance  of  the  filter  is  demonstrated 
on  Figures  26  and  27  which  show  time  histories  of  the  percent 
error  in  X3  and  Xi^ .   As  can  be  seen  from  these  figures,  the 
initial  estimates  of  X3  and  Xi^  are  25  percent  in  error  and 
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these  estimates  are  corrected  to  within  approximately  1  per- 
cent and  3  percent  of  their  respective  true  values.   The 
noise  variances  used  are  those  given  by  the  simulation 
which  generated  the  pitch  angle  data.   The  initial  para- 
meter variances  are  calculated  from  Equation  85. 
Effects  of  variations  in  the  initial  parameter  variances 

The  sensitivity  of  the  filter  to  variations  in  the 
initial  parameter  variances  is  given  in  Table  VIII  and  on 
Figures  28  through  32.   Table  VIII  is  a  summary  of  the 
percent  error  in  the  parameters  of  interest  and  their 
normalized  uncertainties  for  variations  in  the  initial 
parameter  variances  of  -25  percent,  0.0  percent  and  25 
percent  relative  to  the  sample  variances. 

Figures  28  and  29  show  the  variation  of  the  near 
steady-state  errors  in  the  parameters  and  their  normalized 
uncertainties  as  functions  of  the  percent  variations  in 
the  initial  parameter  variances.   Figure  30  shows  the  two 
time  histories  of  the  error  in  the  estimate  of  Xi+  that 
result  for  tv.o  different  initial  values  of  the  parameter 
variances.   v)ne  of  the  initial  variances  is  the  sample 
variance  based  on  the  initial  and  true  values  of  the 
parameters  [Equation  (85)]  and  the  other  is  25  percent 
less  than  this  value. 

The  X3  error  time  histories  are  essentially  identical 
for  both  cases  and  are  not  presented.   The  variations  of 
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Table  VIII.   EFFECTS  OF  VARIATIONS  IN  THE  INITIAL 
PARAMETER  VARIANCES  ON  NEAR  STEADY- 
STATE  PARAMETERS  AND  THEIR  VARIANCES 
(LINEAR  SYSTEM,  MEASUREMENT  ERRORS 
AND  SYSTEM  NOISE) 


Variation  in 
Initial  Parameter 
Variance* 
(Percent) 


Errors  in  Near 
Steady-State 
Parameters 
(Percent) 

Normalized  Near 
Steady-State 
Variances 
(Percent) 

X3        X^ 

^P33        /P^^ 

n         n 

-25           0.72  -0.47  2.42  14.39 

0           0.88  2.66  2.42  15.18 

25           0.98  4.95  2,42  15.81 

relative  to  sample  variance 
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the  parameter  variances  with  time,  for  the  three  initial 
values  considered,  are  shown  on  Figures  31  and  32. 

From  the  information  given  on  the  figures  mentioned 
above,  it  is  obvious  that  neither  X3  nor  P33  are  affected 
by  variations  in  the  initial  variances  within  the  ±25 
percent  range  considered.   The  estimates  of  the  damping 
parameter  as  well  as  its  variance  are  affected,  however, 
by  these  variations.   The  near  steady-state  damping  para- 
meter estimates  vary  almost  linearly  with  initial  variance 
values  from  a  low  of  0.47  percent  to  a  high  of  4.95  percent 
(see  Figure  28)  .   The  variation  in  the  normalized  uncertainty 
is  less  pronounced,  ranging  from  a  low  of  14.39  percent  to  a 
high  of  15. SI  percent  (see  Figure  29).   Nevertheless,  the 
near  steady-state  uncertainties  are  of  sufficient  magnitude 
so  that  the  true  values  of  the  coefficients  are  within 
less  than  one  standard  deviation  of  their  estimated  values 
regardless  of  the  initial  parameter  variance,  when  the 
standard  deviation  is  taken  as  the  square  root  of  the  near 
steady-state  parameter  variance.   Attention  is  also  called 
to  the  fact  that  since  system  noise  is  now  present  in  the 
data,  the  near  steady-state  parameter  variances  are  of 
greater  magnitude  than  was  the  case  when  only  measurement 
errors  were  present  (see  Figures  24,  25,  31  and  32). 
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Filter  response  to  large  errors  in  the  initial  parameter 
estimates 

The  results  of  using  the  extended  Kalman  filter  to 
identify  the  damping  parameter,  Xi+,  when  the  initial  esti- 
mate of  the  parameter  is  in  error  by  100  percent  are  shown 
on  Figure  33  for  two  initial  parameter  variance  estimates. 
When  the  initial  parameter  variance  is  based  on  the  initial 
estimate  of  Xi+ ,  the  near  steady-state  parameter  estimate  is 
approximately  14  percent  high;  for  an  initial  parameter 
variance  that  is  greater  by  approximately  25  percent, 
the  resulting  near  steady-state  value  of  X^  is  nearly  19 
percent  high.   The  corresponding  parameter  variances  1.5 
seconds  after  the  filter  begins  integrating  are  such  that 
the  estimated  values  of  X^  are  within  approximately  one 
standard  deviation  of  the  true  value. 
Comparison  of  filter  performance  to  Chapman-Kirk  technique 

The  parameters  of  interest  and  their  uncertainties 
were  extracted  from  the  given  noisy  data  set  using  the 
previously  described  technique  of  Chapman  and  Kirk  in  order 
to  be  compared  to  the  filter  results.   The  percent  error 
in  the  parameters  and  their  normalized  estimated  standard 
deviations  obtained  with  the  Chapman-Kirk  program  are 

X3:12.5  percent  (1.8  percent) 
Xl^:    3.1  percent  (0.3  percent) 
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The  above  results  are  similar  to  those  presented  in  Chapter 
III  (Table  VI)  in  that  the  estimated  parameter  uncertain- 
ties are  not  of  sufficient  magnitude  to  reflect  the  true 
error  in  the  extracted  parameters.   Whereas  the  extended 
filter  estimates  the  parameters  and  variances  so  that  the 
true  value  of  the  parameter  is  generally  within  one  standard 
deviation  of  the  estimate,  it  is  shown  that  such  is  not  the 
case  when  using  the  Chapman-Kirk  technique.   The  errors  in 
the  above  parameters  are  also  slightly  greater  than  those 
obtained  with  the  extended  filter,  although  this  is  probably 
of  less  consequence  than  the  differences  in  the  uncertainties 
computed  by  the  two  techniques. 

Nonlinear  System  with  Measurement  Errors  Only 

The  pattern  used  in  the  previous  two  sections  of 
generating  a  data  set,  investigating  the  basic  response 
of  the  filter  in  identifying  the  parameters  of  interest 
and  checking  the  sensitivity  of  the  filter  to  variations 
in  the  initial  parameter  variances  is  essentially  repeated 
here  for  the  more  complicated  system  model 

d  +  (C++C6a2)a  +  CC3+C5a2)a  =  0         .  (88) 

In  addition,  ;>ome  effects  of  errors  in  the  measurement 
error  variance  are  also  included  and  discussed.   Some 

comparisons  af  results  achieved  with  tlie  extended  filter 

and  the  Chapav:in-Kirk  technique  are  also  made. 
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The  standard  deviation  of  the  measurement  errors  in 
the  generated  data  is  0.00569  radian.   The  mean  value  of 
the  errors  is  zero. 
Basic  filter  performance 

The  tinite  histories  of  the  four  parameters  of  interest 
as  computed  by  the  extended  filter  while  analyzing  the  data 
just  described  are  shown  on  Figures  34,  35,  36  and  37. 
The  initial  errors  in  the  parameters  are  -25  percent  for 
X3,  Xi^    and  X£  and  -5  percent  for  X5.   The  measurement 
error  variance  is  that  obtained  from  the  simulation.   The 
initial  values  of  the  parameter  variances  are  based  on 
the  difference  between  the  true  values  of  the  parameters 
and  the  above  mentioned  initial  estimates. 

The  reason  the  initial  estimate  of  X5  is  only  5  percent 
low  instead  of  25  percent,  as  is  the  case  with  the  other 
parameters,  is  due  to  its  magnitude  and  a  weakness  in  the 
numerical  scheme.   The  true  value  of  X5  is  1960  based  on 
the  data  in  Table  III  of  Chapter  III.   A  25  percent  initial 
error  in  this  parameter  results  in  an  initial  parameter 
sample  variance  of 

PssCO)  =  (1960- (0.25)1960)2  =  240,000 

which  is  of  sufficient  magnitude  to  cause  the  filter  to 
diverge  for  the  numerical  integration  step  size  of  0.005 
second.   The  divergence  occurs  because  the  large  initial 
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value  of  P55  results  in  large  values  of  P25  and  P15  which 
in  turn  cause  the  values  of  P55  that  are  fed  into  the 
numerical  integration  subroutine  to  be  very  large  [Equa- 
tions (81)  and  (82)].   The  result  of  this  cascading  effect 
is  a  rapid  dl-ecrease  in  P55  in  large  increments  until  it 
becomes  negative  (a  physical  impossibility  since  this  term 
is  the  variamce  of  a  parameter) .   The  structure  of  the 
covariance  equations  is  such  that  a  negative  variance 
causes  numerical  divergence  in  the  solutions  to  several 
of  the  equations. 

To  circumvent  the  problem,  an  initial  error  of  -10 
percent  in  X5,  with  a  corresponding  adjustment  in  P55(0), 
was  tried.   The  result  was  again  divergence.   Finally  an 
initial  esti-siate  of  -5  percent  was  found  to  be  successful. 

The  results  shown  on  Figures  34,  35  and  36  indicate 
that  the  filter  does  an  excellent  job  in  correctly  identi- 
fying the  linear  and  nonlinear  static  restoring  moment 
parameters,  X3  and  X5,  as  well  as  the  linear  pitch  damping 
term,  Xi^.   Tlie  results  for  the  nonlinear  damping  parameter, 
Xg,  are  not  as  gratifying  as  can  be  seen  by  referring  to 
Figure  37.   The  error  in  this  term  is  still  approximately 
-10  percent  after  1.5  seconds.   The  reason  for  this  problem 
is  related  to  the  magnitude  of  this  parameter  as  compared 
to  others  in   the  equation  of  motion.   The  true  value  of 
X3  and  X5  are  160  and  1960  respectively,  although  the 


76 

magnitude  of  X5  is  effectively  reduced  by  at  least  an  order 
of  magnitude  v\fhen  it  is  multiplied  by  a^.   Xt^  and  Xg,  on 
the  other  hand,  have  true  values  of  1.60  and  4.35  respec- 
tively.  After  Xg  is  multiplied  by  a^a  it  is  effectively 
the  smallest  parameter  in  the  equation  of  motion  by  an  order 
of  magnitude  and  it  is  naturally  more  difficult  to  identify, 
since  it  has  the  least  effect  on  the  trajectory. 
Effects  of  variations  in  the  initial  parameter  variances 

The  effects  of  variations  in  the  initial  parameter 
variances  are  summarized  in  Table  IX  and  on  Figures  38 
and  39.   As  can  be  seen  from  these  results,  the  near 
steady-state  values  of  the  parameters  and  their  variances 
are  relatively  insensitive  over  the  range  of  initial  para- 
meter variances  considered.   Only  the  Xg  parameter  estimate 
and  its  variance  show  a  variation  of  more  than  1  percent. 
As  is  obvioiis  from  the  referenced  table  and  figures,  no 
data  are  available  for  initial  parameter  variances  of  50 
percent  above  the  sample  variances;  the  reason  for  this  is, 
again,  that  the  magnitude  of  the  P55  variance  caused  the 
filter  to  diverge. 

AttentiOiT  is  called  to  the  fact  that  for  the  range 
of  initial  variances  considered,  the  near  steady-state 
estimates  of  the  X3,  Xi+  and  X5  parameters  are  all  within 
one  standard  deviation  of  their  true  values  when  the 
standard  deviation  is  taken  as  the  square  root  of  the  near 
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steady-state  variance  for  the  parameter  of  interest.   The 

Xg  parameter  is  always  within  two  standard  deviations  of 

its  true  value. 

Effects  of  errors  in  the  estimate  of  the  measurement  error 

variance 

Up  to  this  point,  exact  values  of  the  measurement 
error  variance,  r,  have  been  used  in  all  runs  made  with  the 
extended  filter.   Exact  knowledge  of  this  quantity  is  avail- 
able because  the  pitch-angle  data  being  used  with  the  extended 
filter  are  computer  generated  with  specified  noise  statistics. 
In  actual  test  situations,  true  values  of  the  noise  para- 
meters may  not  be  known  exactly  and,  in  fact,  methods  for 
determining  both  measurement  error  variances  and  system 
noise  variances  from  noisy  dynamic  data  are  of  current 
research  interest. 

Some  effects  of  errors  in  the  estimates  of  the  measure- 
ment error  variance  on  the  near  steady-state  parameter 
estimates  and  their  uncertainties  are  given  in  Table  X 
and  on  Figures  40  and  41.   These  results  show  that  errors 
in  the  measurement  noise  variance  of  ±25  percent  have 
essentially  no  effect  on  estimates  of  the  four  parameters 
of  interest.   The  uncertainties  of  the  X3,  X^  and  X5 
parameters  are  changed  by  less  than  1  percent  over  the 
range  considered.   Tlie  normalized  Xg  uncertainty  varies  from 
approximately  5  percent  to  12  percent. 
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Filter  response  to  larger  error  in  the  initial  estimate 
of  X5 

The  difficulty  associated  with  numerically  large  initial 
parameter  variance  estimates  has  been  discussed  earlier. 
The  sample  variance  of  P55  which  corresponds  to  a  -5  percent 
error  in  the  initial  X5  parameter  estimate  appears  to  be  an 
approximate  upper  limit  for  the  integration  step  size  and 
physical  constants  used  in  this  analysis.   A  run  was  made, 
hov;ever,  with  the  previously  described  nonlinear  data  set 
containing  only  measurement  errors  where  all  initial  parameter 
estimates  were  in  error  by  -25  percent;  the  P55  variance, 
however,  was  set  equal  to  the  previously  used  sample 
variance  which  corresponds  to  a  -5  percent  initial  estimate 
of  X5.   T?ie  parameter  estimate  time  histories  resulting  from 
tliis  initialization  are  shown  on  Figures  42  through  45. 
These  results  indicate  that  the  near  steady-state  values 
of  the  paraaeters  are  very  nearly  equal  to  those  presented 
when  the  sasTiple  variances  based  on  the  initial  parameter 
estimates  were  used  initially. 
Com.parison  of  filter  performance  to  Chapman-Kirk  technique 

The  linear  and  nonlinear  static  restoring  moment  and 
damping  parameters  were  extracted  from  the  noisy  data  set 
being  considered  in  this  section  using  the  Chapman-Kirk 
program.   TJje  percent  errors  in  the  parameters  and  their 
normalized  estimated  uncertainties  are  given  below: 
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X3:0.36  percent  (0.35  percent) 

Xt,:0.69  percent  (0.43  percent) 

X5:0.41  percent  (0.49  percent) 

X5:4.60  percent  (11.27  percent) 

Comparing  the  above  results  with  those  previously  given 
in  Table  IX  reveals  that  the  estimates  obtained  using  the 
extended  filter  are  slightly  more  accurate  for  all  but  the 
Xi+  parameter.   Since  the  errors  are,  for  the  most  part, 
less  than  1  percent  the  differences  are  essentially 
insignificant.   The  normalized  uncertainties  obtained 
with  both  techniques  are  equally  good. 

Nonlinear  System  with  Measurement  Errors  and  System  Noise 

In  this  section,  as  before,  the  basic  filter  performance, 
sensitivity  of  the  filter  to  initial  parameter  variance  vari- 
ations and  sensitivity  of  the  filter  to  errors  in  the  esti- 
mates of  the  noise  variances  are  investigated.   The  data 
used  in  the  analysis  were  generated  by  numerically  inte- 
grating the  following  nonlinear  equation  with  a  random 
forcing  function, 

a  +  (C4+C5a2)a  +  (C3+C5a2)a  =  w(t)     .         (89) 

Approximately  the  same  noise  variances  used  previously  are 
repeated  once  again;  tlie  variance  of  the  measurement  errors 
is  0.00576  radian  and  the  variance  of  the  system  noise  is 
4.75  rad/sec^. 


82 

Basic  filter  performance 

The  time  histories  of  the  errors  in  the  parameters 
of  interest  obtained  with  the  extended  filter  are  given  on 
Figures  46  through  49.   The  response  time  required  for  the 
parameters  to  reach  their  near  steady-state  values  is 
slightly  greater  than  in  the  simpler  cases  that  have  been 
considered  previously.   The  accuracy  of  the  parameters  is 
similar  to  that  which  was  achieved  for  the  nonlinear  system 
with  only  measurement  errors  in  the  data  with  the  exception 
of  Xi+,  which  is  approximately  5  percent  higher. 
Effects  of  variations  in  the  initial  parameter  variances 

The  effects  of  variations  in  the  initial  parameter 
variances  are  given  in  Table  XI  and  are  shown  on  Figures 
50  and  51.   Only  the  nonlinear  damping  parameter  and  its 
uncertainty  vary  by  more  than  1  percent  for  initial  variance 
variations  of  up  to  50  percent  relative  to  the  sample 
variances. 
Effects  of  errors  in  the  estimates  of  the  noise  variances 

The  effects  of  errors  in  the  estimates  of  r,  the  measure- 
ment error  variance,  and  q,  the  system  noise  variance,  are 
given  in  Table  XII  and  are  shown  on  Figures  52  and  53. 
These  figures  show  the  variation  of  the  near  steady-state 
parameter  errors  and  normalized  parameter  uncertainties  as 
functions  o£  the  error  in  the  estimate  of  the  measurement 
error  variance  for  various  values  of  error  in  the  system 
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noise  variance  estimate.   From  these  results  it  is  obvious 
that  only  Xg  is  affected  by  errors  of  +25  percent  in  the 
estimate  of  the  system  noise  variance  or  measurement  error 
variance . 
Comparison  of  filter  performance  to  Chapman-Kirk  technique 

The  percent  errors  in  the  parameters  and  their  normalized 
estimated  standard  deviations  that  are  obtained  using  the 
Chapman-Kirk  program  with  the  data  currently  being  analyzed 
are  given  below: 

X3:0.62  percent  (0,40  percent) 

Xij:2.63  percent  (1.23  percent) 

X5:0.41  percent  (0.61  percent) 

X5:46.3  percent  (14.8  percent) 

Comparing  these  results  with  the  results  in  Table  XI 
indicates  that  both  techniques  yield  essentially  the  same 
accuracy  in  their  determination  of  X3  and  X5.   The  filter 
does  considerably  better  in  estimating  X5  but  is  slightly 
worse  in  its  estimate  of  Xi+. 

The  uncertainties  computed  using  the  filter  are  of 
sufficient  magnitude  so  tliat  the  true  value  of  the  X3, 
X5  and  Xg  parameters  are  within  one  standard  deviation 
of  tlie  filter  estimates;  the  true  value  of  Xi+  is  within 
two  standard  deviations. 
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When  using  the  Chapman-Kirk  technique,  only  X5  is 
within  one  estimated  standard  deviation  of  its  true  value 
X3  and  Xi+  are  within  two  estimated  standard  deviations 
of  their  respective  true  values.   The  extracted  value  of 
Xg  is  not  within  three  estimated  standard  deviations  of 
its  true  value. 
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Figure  38.   Variation  of  Near  Steady-State  Parameter  Error 
with  Initial  Parameter  Variance  (Nonlinear 
System,  Measurement  Errors  Only) 
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Figure  50.   Variation  of  Near  Steady-State  Parameter  Error 
X'lith    Initial  Parameter  Variance  (Nonlinear 
System,  Measurement  Errors  and  System  Noise) 
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Figure  51,   Variation  of  Normalized  Near  Steady-State 

Parameter  Uncertainty  with  Initial  Parameter 
Variance  (Nonlinear  System,  Measurement 
Errors  and  System  Noise) 
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Figure  52.   Variation  of  Near  Steady-State  Parameter 
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Measurement  Error  Variance  (Nonlinear 
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CHAPTER  VI 
CONCLUSION 

An  analysis  of  methods  for  extracting  aerodynamic 
coefficients  from  dynamic  test  data  has  been  conducted. 
The  primary  concern  of  this  analysis  has  been  to  determine 
the  accuracy  with  which  these  coefficients  can  be  determined 
from  dynamic  data  containing  system  noise  and  measurement 
errors.   Both  linear  and  nonlinear  systems  have  been  con- 
sidered with  emphasis  on  the  latter. 

Two  methods  of  extracting  the  coefficients  of  interest 
and  estimating  their  uncertainties  ha/e  been  analyzed.   The 
first  of  these  is  a  deterministic  tecimique  due  to  Chapman 
and  Kirk.     The  second  technique  considered,  based  on  the 
extended  Kalman  filter"^'^^  with  augmented  state  vector, 
is  stochastic  in  nature  and  has  been  developed  and  applied 
in  this  paper. 

Both  techniques  considered  will  extract  aerodynamic 
coefficients  from  noisy  dynamic  data.   The  degree  to  which 
they  are  successful  has  been  presented  in  detail  in  the 
preceding  chapters. 

The  Chapman- Kirk  Technique 

The  presence  of  noise  in  the  dynamic  data  being  used  by 
the  Chapman-KLrk  technique  affects  the  accuracy  with  which 
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this  technique  can  be  used  to  determine  the  coefficients 
of  interest.   When  measurement  errors  only  are  present  in 
the  data,  these  effects  have  been  found  to  be  negligible 
on  linear  and  nonlinear  static  pitching  moment  coefficient 
derivatives;  the  effects  on  the  linear  and  nonlinear  pitch 
damping  coefficients  are  more  pronounced,  especially  for 
the  higher  measurement  errors.   The  error  in  the  linear 
pitch  damping  coefficient  is  as  high  as  25  percent  in  one 
case.   The  greatest  error  in  the  extracted  nonlinear  pitch 
damping  coefficient  is  227  percent.   There  is  an  obvious 
inconsistency  in  the  extracted  values  of  this  nonlinear 
coefficient  since  another  data  set,  which  also  contained 
measurement  errors  with  similar  statistical  properties, 
yielded  a  value  for  this  term  which  was  in  error  by  only 
5  percent. 

The  estimation  of  the  uncertainties  of  the  extracted 
coefficients  vvhen  the  data  contained  only  measurement  errors 
is  quite  adequate.   In  every  case  considered,  the  true  value 
of  the  coefficient  of  interest  was  found  to  be  within  three 
estimated  standard  deviations  of  the  extracted  value.   The 
estimated  standard  deviations  of  a  given  coefficient  are  also 
consistent  for  different  data  sets  with  the  same  measurement 
error  statistics. 

Use  of  data  containing  measurement  errors  and  system 
noise  has  a  more  adverse  effect  on  the  accuracy  of  coeffici- 
ents determined  with  tlie  Chapman-Kirk  technique.   As  is  the 
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case  when  analyzing  data  containing  only  measurement  errors, 
the  pitch  da:mping  coefficients  are  the  ones  most  affected. 
The  errors  in  these  coefficients  range  from  1  percent  to 
43  percent  for  the  linear  term  and  45  percent  to  600  percent 
for  the  nonlinear  term  over  the  range  of  system  noise  con- 
sidered. 

Perhaps  of  even  greater  consequence  than  the  errors  in 
the  damping  parameters  is  the  fact  that  their  estimated 
uncertainties  are  too  small  to  indicate  the  true  errors  in 
the  coefficients  when  the  data  contain  both  system  noise  and 
measurement  errors.   This  is  a  direct  consequence  of  using 
the  root-mean-square-error  between  the  final  solution  to 
the  equation  of  motion  and  the  data  to  calculate  the  parameter 
uncertainties.   This  error  term  is  an  adequate  representation 
of  noise  in  the  data  only  if  all  errors  are  randomly  super- 
imposed on  the  pitch  angle  output  after  the  equation  of 
motion  has  been  integrated.   The  integrated  effects  of  the 
random  accelerations  representing  system  noise,  however, 
are  such  that  the  final  errors  appearing  in  the  pitch  angle 
due  to  the  system  noise  are  not  necessarily  random,  or 
independent  from  one  discrete  time  to  another. 

Extended  Kalman  Filter 

The  feasibility  of  using  an  extended  Kalman  filter  with 
parameter  augmented  state  vector  for  determining  the  values 
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o£  aerodynamic  coefficients  and  their  uncertainties  from 
dynamic  test  data  has  been  demonstrated  for  a  one-degree- 
of-freedom  system. 

The  specific  filter  used  here  generally  reaches  near 
steady-state  estimates  of  the  parameters  in  less  than  one 
second  for  the  system  model  and  error  combinations  considered. 

For  linear  systems  with  measurement  errors  only  in  the 
data,  the  extended  filter  yielded  estimates  of  both  the 
static  pitc5iing  moment  coefficient  derivative  and  the  pitch 
damping  coefficient  to  within  1  percent  of  their  respective 
true  values.   Slightly  less  accurate  determinations  of  the 
same  two  parameters  are  obtained  using  the  filter  when 
both  measuremt^nt  errors  and  system  noise  are  present  in  the 
data.   Variations  in  the  initial  parameter  variances  of 
25  percent  seem  to  affect  essentially  only  the  pitch  damping 
term  and  then  only  when  both  measurement  errors  and  system 
noise  are  present  in  the  data. 

For  nonlinear  systems  with  measurement  errors  only  in 
the  data  the  estimates  obtained  with  the  extended  filter 
for  the  linear  and  nonlinear  static  pitching  moment  co- 
efficient derivatives  as  well  as  the  linear  pitch  damping 
coefficient  are  within  approximately  1  percent  or  less 
of  their  respective  true  values  for  a  variety  of  initial 
parameter  variances  and  noise  variance  estimates.   The 
error  in  the  estimate  of  the  nonlinear  damping  coefficient 
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varies  between  11  percent  and  14  percent.   None  of  the 
coefficients  are  particularly  sensitive  to  variations  in 
the  initial  jaarameter  variances  or  noise  variances  con- 
sidered.  As  is  the  case  for  the  linear  system,  the 
extended  flitter  is  less  accurate  in  its  parameter  identi- 
fication whe-Ta  both  measurement  errors  and  system  noise 
of  the  magnitudes  considered  here  are  present  in  the  data. 
The  linear  d'^mping  coefficient  is  most  affected  by 
having  the  aiddition  of  system  noise  in  the  data.   The 
error  in  this  coefficient  increases  from  approximately 
1  percent  to  approximately  6  percent  for  the  cases  con- 
sidered.  Erxors  in  the  noise  variance  estimates  of 
+  25  percent  2a,ave  little  or  no  appreciable  effect  on  the 
linear  and  nonlinear  static  pitching  moment  coefficient 
derivatives  or  the  linear  pitch  damping  coefficient,  and 
only  a  slighl  effect  on  the  nonlinear  pitch  damping 
coefficient. 

For  the  range  of  initial  parameter  errors,  initial 
parameter  variances  and  noise  variances  considered,  the 
extended  KaLasian  filter  produces  an  excellent  estimate  of 
the  parameter  uncertainties  after  integrating  the  Ricatti 
equations  fcx  approximately  1.5  to  2.0  seconds.   For  every 
case  considcsared  where  the  initial  parameter  estimates  are 
within  25  percent  of  their  true  values,  the  near  steady- 
state  parameter  value  is  within  two  standard  deviations 
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of  its  true  A-alue,  and  in  the  majority  of  cases  within 
one  standard  deivation,  when  the  standard  deviation  is 
taken  to  be  the  square  root  of  the  near  steady-state 
parameter  variance. 

Comparison  of  the  Two  Techniques 

There  is  a  natural  tendency  when  analyzing  two  tech- 
niques to  compare  the  results  achieved  when  solving  similar 
problems.   Both  techniques  investigated  have  strengths  and 
weaknesses. 

From  the  preceding  chapters  it  is  obvious  that  the 
Chapman-Kirk  technique  requires  less  information  initially 
than  the  extended  Kalman  filter  since  it  is  only  necessary 
to  make  initial  estimates  of  the  parameters  in  the  former 
whereas  the  parameters,  their  variances  and  the  noise 
variances  must  all  be  estimated  initially  when  using  the 
extended  filter. 

There  is  no  provision  to  weigh  the  initial  parameter 
estimates  when  using  the  Chapman-Kirk  technique  whereas 
this  is  done  through  the  initial  variances  when  using 
the  extended  filter.   This  feature  is  particularly 
important  if  certain  parameters  are  known  with  signifi- 
cantly higher  accuracy  before  the  extraction  process  begins. 
The  Chapman-Kirk  technique  will  sometimes  actually  down- 
grade good  initial  estimates  when  analyzing  noisy  data 
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whereas  the  extended  filter  will  tend  to  improve  the 
estimate  or  at  least  not  downgrade  it  for  the  initial 
parameter  variances  considered. 

From  an  economic  point  of  view^,  both  techniques  are 
essentially  equal  concerning  computational  costs.   The 
extended  filter  must  numerically  integrate  more  equations 
bu  only  has  to  do  this  once;  the  Chapman-Kirk  technique, 
being  an  iterative  scheme,  must  integrate  fewer  equations 
more  times. 

The  programming  requirements  of  the  extended  filter 
are  less  than  those  of  the  Chapman-Kirk  technique  since 
the  former  requires  no  matrix  inversions,  simultaneous 
equation  solutions  or  logic  for  including  or  excluding 
certain  parameters  from  the  extraction  process. 

Final  Remarks 

The  research  reported  in  this  dissertation  is  part 
of  a  larger  and  continuing  research  program  whose  goal 
is  a  critical  study  of  methods  for  extracting  aerodynamic 
coefficients  from  dynamic  data.   Since  the  research  is 
of  a  continuing  nature  it  seems  appropriate  to  mention 
some  areas  where  additional  work  should  be  contemplated. 

The  results  presented  here  are  entirely  for  a  one- 
degree-of- freedom  system.   The  effect  of  system,  noise  and 
measurement  errors  on  coefficients  extracted  from  coupled, 
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multiple  degree  of  freedom  systems  may  prove  to  be 
extremely  interesting  and  is  the  next  logical  area  to  be 
investigated. 

The  response  of  the  extended  filter  to  highly  erroneous 
initial  parameter  estimates  seems  quite  good  for  the  limited 
nuniber  of  cases  investigated.   It  seems  possible  to  take 
advantage  of  this  feature  by  constructing  an  extraction 
algorithm  which  employs  both  an  extended  filter  and  a 
Cht.pman-Kirk  or  similar  technique  in  tandem  with  the 
extended  filter  providing   he  required  accurate  initial 
estimates  for  the  differential  corrections,  iterative 
technique . 

Additioi:al  work  in  relating  the  effects  of  inadequate 
or  erroneous  "Modeling  to  the  noise  problem  should  also 
prove  to  be  of  interest. 

Finally,  when  considering  new  techniques  such  as  the 
use  of  the  extended  Kalman  filter,  it  certainly  must  be 
concluded  that  additional  numerical  investigations  will 
certainly  add  to  the  confidence  or  more  clearly  define 
the  limitations  of  this  method  of  coefficient  extraction. 


APPENDIX  I 

CHAPMAN-KIRK  COEFFICIENT  EXTRACTION  NOMENCLATURE 
LIST  AND  PROGRAM  LISTING 


Nomenclature  List 


PROGRAM 
VARIABLE 

H 

I  TO 

TMAX 

TZERO 

N 

X2(l) 
XZ(2) 

xzrs) 

XZ(4) 

XZ(5) 


MATH  SYMBOL 


max 


N 


ao 


«0 


3a 


3C- 


d   9a 
dt  9  0  1 


3a 


3C2 


t  =  0 


t=0 


t  =  0 


DEFINITION 

Numerical  integration 
step  size  (sec) 

Frequency  of  numerical 
integration  output 

Cutoff  time  for  numeri- 
cal integration  (sec) 

Initial  time  for 
numerical  integration 
(sec) 

Number  of  first  order 
differential  equations 
to  be  integrated 

Initial  condition  of  a 
(rad) 

Initial  condition  of  d 
(rad/sec) 

Initial  condition  for 
parametric  differential 
equation 

Initial  condition  for 
parametric  differential 
equation 

Initial  condition  for 
parametric  differential 
equation 
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PROGRAM 
VARIABLE 

XZ(6) 


XZ(7) 


MATH  SYMBOL 


XZ(8) 


XZ(9) 


XZ(IO) 


XZCll) 


XZ(12) 


XZ(15) 


XZ(14) 


XZ(15) 


XZ(16) 


d   3a 


dt  8C2 


3a 


3C- 


d   3a 
^  9C3 


9a 


3Ci 


d   3a 


dt  3C4 


3a 


3C. 


d   3a 
^  3C5 


3a 


3C( 


d   3a 
3T  3C6 


t=0 


t=0 


t  =  0 


t  =  0 


t  =  0 


t=0 


t  =  0 


t  =  0 


t  =  0 


3_a_ 

3C7 


d   3_a_ 
dT  3C7 


t  =  0 


DEFINITION 

Initial  condition  for 
parametric  differential 
equation 

Initial  condition  for 
parametric  differential 
equation 

Initial  condition  for 
parametric  differential 
equation 

Initial  condition  for 
parametric  differential 
equation 

Initial  condition  for 
parametric  differential 
equation 

Initial  condition  for 
parametric  differential 
equation 

Initial  condition  for 
parametric  differential 
equation 

Initial  condition  for 
parametric  differential 
equation 

Initial  condition  for 
parametric  differential 
equation 

Initial  condition  for 
parametric  differential 
equation 

Initial  condition  for 
parametric  differential 
equation 
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PROGRAM 
VARIABLE 

XZ(17) 


XZ(18) 

Q 

V 

AMI 


MATH  SYMBOL 


3a 


A 
D 

CM/^O 

CMA2 

CMA4 

CMQO 
CMQ2 
CMDA 


3C< 


d  da 

3T  dCg 


A 
d 

■'mao 
■'ma  2 

"ma  4 

"mqo 
"mq2 

'mSa 


t  =  0 


t=0 


DEFINITION 

Initial  condition  for 
parametric  differential 
equation 

Initial  condition  for 
parametric  differential 
equation 

Dynamic  pressure 
(lb/ft2) 

Freestream  velocity 
(ft/sec) 

Moment  of  inertia  about 
an  axis  through  the 
vehicle  C.G.  and  normal 
to  the  pitch  plane 
(slug  .  ft2) 

Vehicle  reference  area 

(ft2) 

Vehicle  reference 
length  (ft) 

Static  pitching  moment 
coefficient  derivative 
(rad-i) 

Static  pitching  moment 
coefficient  derivative 
(rad-3) 

Static  pitching  moment 
coefficient  derivative 
(rad-5) 

Pitch  damping  coeffi- 
cient 

Pitch  damping  coeffi- 
cient (rad~2) 

Static  pitching  moment 
coefficient  at  a  =  0 
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PROGRAM 
VARIABLE 

M 

JJJMAX 


MATH  SYMBOL 


M 


DEFINITION 

Number  of  experimental 
data  points 

Maximum  number  of 
iterations  allowed 
before  program  termi- 
nation 


EP 

AE(I) 

CI 
C2 

C3 


a    (tj 
exp  ^  ' ^ 


Ci 


Convergence  criteria 
for  change  in  RMS 
ERROR  (rad) 

Experimental  values  of 
a  (rad) 


C    q  Ad 

mao^ 


C4 


C    q  Ad2 
mqo  ^ 

2VI 


C5 


C6 


C^  9^  Ad 

I 


C^n9^  Ad2 

mq2  ^ 

2VI 


C7 


C7 


C    q  Ad 


C8 


C^^  q  Ad 
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Categories  of  Input  Data 

Numerical  integration  constants .- -These  constants 
include  the  numerical  integration  step  size,  frequency  of 
output,  time  at  which  the  numerical  integration  is  to  stop, 
initial  time,  and  the  number  of  first  order  equations  to 
be  integrated. 

Initial  condition^ .- -These  data  are  the  first  estimates 
of  the  initial  conditions  for  the  equation  of  motion  as  well 
as  the  initial  conditions  for  the  parametric  differential 
equations . 

Aerodynamic  and  physical  constants .- -These  constants 
are  the  dynamic  pressure  and  freestream  velocity  that 
existed  when  the  experimental  pitch  angle  data  were  re- 
corded, the  moment  of  inertia  about  an  axis  normal  to  the 
pitch  plane  and  through  the  center  of  gravity  of  the  vehicle 
of  interest,  a  vehicle  reference  area,  and  a  vehicle  refer- 
ence length. 

Aerodynamic  coefficient  estimates .- -These  data  are 
the  first  estimates  of  the  aerodynamic  coefficients  to  be 
extracted. 

Convergence  criteria .- -These  constants  are  the  number 
of  experimental  data  points  to  be  read,  the  maximum  number 
of  iterations  to  be  allowed,  and  the  tolerance  against 
which  the  root-mean-square  error  difference  is  tested  for 
convergence. 
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Experimental  data. --These  data  are  the  experimental 
pitch  angle  v'alues  from  which  the  aerodynamic  coefficients 
are  to  be  extracted.   The  program  assumes  the  time  between 
each  data  point  is  constant,  and  this  time  must  be  equal 
to  the  product  of  the  numerical  integration  step  size  and 
frequency  of  numerical  integration  output  that  were  given 
on  the  first  data  card. 
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APPENDIX  II 

EXTENDED  KALMAN  FILTER  NOMENCLATURE  LIST 
AND  PROGRAM  LISTING 


Nomenclature  List 


PROGRAM 
VARIABLE 

H 

TZERO 

TMAX 

ITO 

N 

M 
XZ(I) 

Q 

R 

ZCD 
x(i) 


MATH  SYMBOL 


At 


'max 


N 


M 


DEFINITION 

Numerical  integration 
step  size  (sec) 

Initial  time  for  num- 
erical integration  (sec) 

Cutoff  time  for  numerical 
integration  (sec) 

Frequency  of  numerical 
integration  output 

Number  of  first  order 
equations  to  be  integrated 

Number  of  data  points  • 

Initial  conditions  for 
state  variable  and 
variance  estimates 

System  noise  variance 
(rad^/sec"^) 

Measurement  error 
variance  (rad^) 

Experimental  data 
values  of  a 

Estimated  value  of 
a  (rad) 
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PROGRAM 
VARIABLE 

XC2) 


X(3) 

X(4) 

X(5) 

X(6) 

X(7) 

X(8) 

X(9) 

X(10) 

X(ll) 

X(12) 

X(13) 

X(14) 

X(15) 

X(16) 

X(17) 

X(18) 

X(19) 

X(20) 

X(21) 

X(22) 

XC23) 

X(24) 


MATH  SYMBOL 


1  1 


12 


1  3 


1  h 


1  5 


16 


22 


23 


2k 


25 


26 


3  3 


31+ 


35 


36 


4  1+ 


4  5 


kS 


DEFINITION 

Estimated  value  of  a 
(rad/sec) 

Estimated  value  of  C3 

Estimated  value  of  C^ 

Estimated  value  of  C5 

Estimated  value  of  Cg 

Variance  of  X^ 

Covariance  of  X^  and  X. 

Covariance  of  X^  and  X. 

Covariance  of  X^  and  Xj 

Covariance  of  X^  and  X, 

Covariance  of  X^  and  X, 

Variance  of  Xo 


Covariance  of  X2  and  X3 

Covariance  of  X2  and  X^ 

Covariance  of  X2  and  X5 

Covariance  of  X,  and  X^ 
Variance  of  X3 

Covariance  of  X,  and  X^^ 

Covariance  of  X3  and  X5 

Covariance  of  X3  and  Xg 
Variance  of  Xr 


Covariance  of  X  ^^  and  X^ 
Covariance  of  X^  and  Xg 
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PROGRAM 

VARIABLE 

MATH  SYMBOL 

X(25) 

P55 

X(26) 

P56 

X(27) 

Pse 

A(I) 

• 

X 

DEFINITION 


Variance  of  X^ 

Covariance  of  X^.  and  X^ 

Variance  of  X, 

Derivative  of  X  with 
respect  to  time 


145 
Categories  of  Input  Data 

Numerical  integration  constants .- -These  constants  are 
the  numerical  integration  step  size,  frequency  of  output, 
time  at  which  the  integration  is  to  stop,  initial  time,  the 
number  of  first  order  equations  being  integrated  and  the 
number  of  data  points  to  be  read. 

Initial  conditions .- -These  data  are  the  initial  estimates 
of  the  state  variable  components,  state  variable  variances  and 
ccvariances . 

Noise  variance  estimates .- -These  two  entries  are 
the  estimates  of  the  measurement  error  variance  and  system 
noise  variance. 

Pitch  angle  data. --These  data  are  the  experimental 
values  of  the  pitch  angle  from  which  the  filter  estimates 
the  states  and  parameters.   The  time  between  data  points 
is  not  input  but  is  assumed  to  be  tlie  same  as  the  numerical 
integration  step  size  input  on  the  first  data  card. 
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